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1.  INTRODUCTION 


In  recent  work,  referred  to  as  I’Zl  (Powell  and  Zielinski  1989),  we  developed  a  steady- 
stale  model  for  investigating  the  properties  of  plasma  discharges  in  an  cleclrolheniial  tl'/l  ) 
gill).  A  numln'r  of  experimental  measurements  were  also  discussed  and  a  comparison  of  the 
ex|>eiimenlnl  and  theoretical  results  was  undertaken.  In  this  report,  we  will  be  interested 
primarily  in  extending  the  calculations  to  the  time-dependent  regime.  We  will  also  account 
for  certain  nonidcal  eflects,  neglected  in  our  previous  work,  but  which  can  be  important 
for  plasmas  in  the  temperature  and  density  range  characteristic  of  these  discharges.  More 
<;xteiisiv<’  experimental  me;isurenicnts  will  also  he  presented  and  comimred  with  the  calcu¬ 
lations.  It  is,  in  fact,  these  additional  measurements  and  their  lack  of  agreement  with,  the 
steady-stale  model  that  have  motivated  the  lime-dependent  calculations. 

1  he  intent  of  both  the  theoretical  and  experimental  work  described  here  is  to  provide 
some  insight  into  the  nature  of  the  a-pillary  discharge.  It  is  hoped  that  ultimately  a  model 
can  he  developed  which  can  be  couijlcd  to  interior-ballistic  codes,  as  well  as  to  models  for  the 
appropriate  power  supplies,  and  provide  a  complete  description  of  the  electrothermal  gun. 
"I  he  cornjjlexily  of  the  analysis  as  a  wliole  has  been  the  barsic  motivation  for  various  groups 
conccnti ating  on  a  single  aspect  of  the  problem. 

During  the  la.st  decade  or  .so,  considerable  effort  has  been  expended  in  the  study  of 
problems  similar  to  the  capillary  discharge.  This  emphasis  has  arisen  because  of  the  impor¬ 
tance  of  ablation-stabilized  arcs  for  a  variety  of  applications.  Probably  the  most  complete 
e.'<|)eriiiK'ntal  work  was  (he  early  investigations  hy  Ibrahim  (1980)  who  undertook  quite  an 
extensive  set  of  diagnostic  measurements  of  various  arc  |)roperties.  A  fairly  simj)le  theoret¬ 
ical  model  Wins  also  |»ro|)(j.sed  and  used  to  analyze  the  experimental  data.  Later,  Kovitya 
and  bowke  (1!)8'1)  undertook  a  thoroughgoing  theoretical  treatment  of  the  problem.  Their 
model  was  one-dimensional  and  i)redicted  the  pioj)erli<!s  of  the  discharge  under  steady-state 
condilioiiH.  UcHults  of  the  calculation.s  were  comj)ared  to  a  very  limited  extent  with  the 
expel imeiital  wc^rk  of  Ibrahim  (1980),  and  also  of  Niemeyer  (1978),  and  reasonabK’  good 
agreement  between  theory  and  cxj;erim'’nt  wa.s  claimed.  More  recent  work,  both  theoretical 
and  experimental,  ha.s  been  carried  out  by  llucl.ti  and  .Niemeyer  (1980).  With  regard  to  the 
specihe  problem  (>{  electi'olhcrmal-guii  cajnllaries,  models  of  varying  complexity  have  been 
projio.sed  hy  Tidmaii,  'I'liio,  (Johlslein,  and  Spicer  (1980),  hy  Loeb  and  Kaplan  (1989),  by 
(  liiyssomalliH,  Mariiux;,  Ricci,  and  Cook  (1988),  by  Mohanti,  Gilligan,  and  Bourham  (,Mo- 
haiili  and  (lilligan  1990;  Mohanti,  (jilligan,  and  Bourrnan  1991),  and  by  various  investigators 
at  Saiidia  National  Laboialorieh  (Woodliii  1991).  These  treatments  range  from  zero  dimen- 
hiumvl  and  steady  stale,  to  more  sojihisticaled  lime-dependent,  one-dir  icnsional  calculations 
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which  arc  coupled  lo  belli  llie  iiulsc-forniing  network  and  to  inlcrior-ljallisf ic  modch  of  the 
working  fluid. 

d'he  model  and  calculations  both  here  and  in  oiir  prcviou.s  work  arc  similar  to  tin'  oiic- 
diinen.sioiial  ticalments  above.  Our  inlcniion  has  been  to  develop  a  capaliility  in  IfKI.  for 
describing  Lite  ET  capillary  in  order  lo  provide  basic  data  for  modeling  the  interior  balli.stics 
of  these  de\'ices.  Furthermore,  we  have  available  models  and  comimtcr  code.s  that  were 
developed  in  conjunction  with  the  study  of  plasma  arcs  in  railguns.  These  codes  solve  the 
basic  magnetohydrodynamic  conservation  equations  for  the  plasma,  as  well  as  calculate  its 
conductivity,  its  ionization  state,  and  approximate  pertinent  equations  of  state.  This  report 
will  basically  describe  our  first  efforts  at  modifying  these  models  and  programs  lo  treat  the 
ET-cai)illary  problem,  as  well  as  present  some  preliminary  experimcnlal  measurements  of 
the  plasma  jiropcrties. 

The  organization  of  the  report  is  as  follows.  In  Sec.  2,  we  describe  the  model  and  major 
assumptions.  In  Sec.  3,  the  theoretical  formalism  necessary  for  treating  the  problem  is 
devclojied.  In  Sec.  4,  wc  di.scuss  a  simple,  limiting-ca.se  calculation;  the  simplified  model  is 
one  which  is  isothermal  and  steady  state.  In  Sec.  5,  the  experimental  measurements  alluded 
to  above  are  described  in  some  detail.  Section  6  contains  results  of  the  calculations  and  a 
comparison  of  theoretical  and  experimental  results.  Finally,  in  Sec.  7,  we  attempt  lo  discuss 
rather  carefully  the  limitations  of  the  calculations,  likely  reasons  for  the  lack  of  agreement 
between  some  aspects  of  calculations  and  experiment,  and  some  ways  in  which  both  the 
model  and  experiments  can  be  improved  in  future  work. 

2.  MODEL  AND  ASSUMPTIONS 

The  basic  model  that  will  be  cm[)loycd  in  the  calculations  is  shown  schematically  in  Fig.  1. 
A  potential  difference  applied  between  the  anode  and  cathode,  located  at  c  =  0  and  ^  =  f  in 
the  figure,  produces  a  current  i  which  is  conducted  through  a  plasma  arc  in  the  interior  of 
the  cylinder  shown.  The  current  olimically  heats  the  plasma  to  temperatures  of  several  tens 
of  thousands  of  degrees.  Radiation  flux  q  from  the  heated  plasma  then  ablates  additional 
material  from  the  capillary  wall.  This  material  replaces  the  plasma  which  is  shown,  via  the 
curved  arrows,  flowing  out  the  end  of  the  capillary  lube  located  at  c  =  The  left  end  of  the 
tube,  located  at  ^  =  0,  is  closed.  The  radius  of  the  capillary,  not  shown  in  the  figure,  is  r;,. 

In  PZ*,  we  undertook  calculations  with  this  model  in  which  it  was  assumed  that  the 
plasma  could  be  treated  in  the  steady-state  approximation.  In  the  present  work  we  will 
extend  tlio.se  calculations  to  the  time-dependent  regime.  In  addition,  we  will  include  certain 


Figure  1.  Schematic  Diagranj  of  Capillary  Model. 

nonideal  effects  which  can  occur  for  plasmas  in  this  temperature  and  density  range.  These 
nonideal  effects  are  important  whenever  the  density  of  the  plasma  is  sufficiently  high  that, 
because  of  nonncgligible  Coulomb  interactions,  the  potential  energy  is  not  negligible  relati\'e 
to  the  thermal  energy.  Nonideal  plasmas  have  received  considerable  attention  in  the  last  few 
years,  and  there  now  exist  an  enorm.ous  number  of  models  (See,  e.g.,  Gunther  and  Radtke 
1981)  for  calculating  the  effects  of  these  important  interactions. 

Nonideal  effects  can  be  expected  to  produce  the  following  consequences.  First,  and 
probably  most  important,  is  a  reduction  in  the  plasma  conductivity  at  a  fixed  temperature 
and  pressure.  Roughly  speaking,  this  reduction  occurs  because  short-range  collisions  between 
electrons  and  ions,  which  occur  infrequently  and  are  negligible  in  low-density  plasmas,  cannot 
be  neglected  in  the  nonideal  case.  To  correct  the  conductivity,  we  have  adapted  a  heuristic 
method  developed  by  Zollweg  and  Licbermann  (1987),  the  principal  result  of  which  is  to 
replace  the  Debye  shielding  parameter  in  the  us’.  al  Spitzer  expression  (Spitzer  1965)  for 
conductivity  by  an  effective,  larger  value.  The  second  consequence  of  nonncgligible  Coulomb 
interactions  is  that  the  nominally  ’Tree”  electrons  in  the  plasma  are  no  longer  completely 
free,  but  remain  weakly  bound  to  the  ions.  Consequently,  there  is  an  effective  reduction  in 
the  ionization  potential  associated  with  each  species,  since  less  energy  is  required  to  remove 
an  electron  to  a  weakly  bound  state  than  to  a  completely  free  state.  In  addition,  one  can 
obviously  expect  a  contribution  to  the  internal  energy  of  the  pleisma  from  the  Coulomb 
interaction  and  a  resulting  change  in  the  pressure  as  a  function  of  temperature  and  density. 


3 


Since  tfic  Coulomb  interaclion  is  negaiivc.  the  pressure  is  reduced  from  that  wliicli  would 
be  obtained  in  the  ideal  case.  Since  the  change  in  both  the  energy  and  pressure  tend  to 
be  very  small  indeed,  at  least  for  cases  of  interest  to  us,  we  have  neglected  these  particular 
nonideal  effects.  W’e  do,  however,  account  for  the  lowering  of  ionization  potential.  .A  similar 
approximation  has  been  employed  for  railgun  plasma  armatures  by  Rolader  and  Balteh 
(19S9).  and  we  will  use  the  sanu'  model  for  the  ionization  reduction  as  they  have  employed. 
The  model  is  that  developed  by  Ebeling  and  Sandig  (1973). 

Two  other  minor  improvements  have  been  made  in  the  calculations  of  PZl.  First,  we  have 
attempted  to  improve  the  computation  of  the  electronic  partition  functions  of  carbon  and 
its  ions  by  including  additional  terms  in  the  series.  Second,  we  have  included  a  calculation 
of  the  electronic  excitation  energy  and  included  this  contribution  in  the  internal  energy  of 
the  plasma.  Both  of  these  changes  are  di.scussed  in  greater  detail  in  the  ensuing  pages,  but 
neither  makes  a  particularly  significant  change  in  the  results  of  PZl. 

Other  approximations  and  assumptions  that  were  employed  in  PZl  are  also  used  in  this 
calculation:  The  model  is  one  dimensional  so  that  it  is  assumed  that  all  flow  variables 
are  nearly  constant  across  the  cross  section  of  the  capillary,  except  possibly  within  a  small 
distance  of  width  S  near  the  walls.  This  region  is  designated  as  the  ’’boundary"  and  is 
so  labelled  in  Fig.  1.  In  all  final  equations  we  take  the  limit  as  (5  — »  0.  It  is  assumed 
that  the  plasma  consists  of  a  homogeneous  mixture  of  hydrogen  and  carbon  neutrals  and 
ions,  which  resulted  from  dissociated  polyethylene,  the  substance  assumed  to  compose  the 
solid  capillary.  The  carbon  is  assumed  to  be  at  most  doubly  ionized.  We  neglect  both 
turbulent  and  magnetohydrodynamic  effects.  The  latter  approximation  is  reasonable  since 
the  magnetic  pressures  turn  out  to  be  small  relative  to  the  hydrostatic  pressures;  the  neglect 
of  turbulence  is  not  so  easily  justified  and  requires  further  study. 

3.  THEORETICAL  FORMALISM 

3.1  Derivation  of'  Conservation  Equations  with  Ablation.  The  most  straight¬ 
forward  method  for  obtaining  the  one-dimensional  equations  which  express  conservation  of 
mass,  momentum,  and  energy  and  which  account  for  ablation  is  to  average  the  appropriate 
two-dimensional  equations  over  the  cross  section  of  the  capillary.  We  therefore  begin  the 
analysis  with  the  standard  axisymmetric  two-dimensional  results  in  cylindrical  coordinates. 
Let  />,  P,  and  E  be  the  density,  pressure,  and  specific  internal  energy  of  the  plasma.  Let  v 
be  velocity  which  has  r  component  u  and  z  component  u;.  Let  J,  cr,  and  q  be  the  current 
density,  electrical  conductivity,  and  heat  flux.  Then,  expressed  in  conservation  form,  the 
equations  can  be  written 
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(1) 


dp  15,  ^  d  ^ 

*  +  r  a;*’''"'* "  a;*'’"'*  "  “ 

5^.15  ^  '  2^ 


(2) 


and 


Ikc  +  fV2)|  +  ;|:|r/«<(£  +  iV2)l  +  £|;>M£  +  '’V2)|  +  +  ^(Pu,)  = 


15  5o. 

- irirqr)--^.  (3) 

a  T  or  az 


\Vc  now  assume  that  all  the  plasma  variables  arc  nearly  constant  across  the  arc  cross 
section  except  very  near  the  boundary.  We  then  multiply  through  each  of  the  equations 
above  by  2nrdr  and  integrate  from  r  =  0  to  r  =  rj,  —  <5,  i.  e.,  up  to  the  inside  edge  of  the 
boundary  layer.  The  spatial  average  of  any  variable  F  is  defined  as 


2n 

^  =  T  [  rF{r,z)dr  (4) 

A  Jo 

where  A  is  the  cross  sectional  area  given  by  A  =  ffrf.  In  obtaining  the  above  result  we  have 
taken  6  to  be  small  relative  to  Vh-  We  furthermore  assume  that  since  the  variables  are  nearly 
constant  that  the  average  of  all  products  can  be  replaced  by  the  product  of  averages.  The 
above  integrations  then  yield 


and 


5p  5  2 

—  +  - - PhUk  =  Pa 

at  az  Tk 


5  0  2 

TTlPU.’)  +  —[pW  )  +  —{pkUkWk) 

at  az  Tk 


dz 


(r>) 


^[p(£:  +  1)V2)]  +  +  v^/2)  +  ^  “ 

^PkUk(Ek  4  vl/2)  -  —Qrk-  (7) 

r-fc  rk 
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llrrc  the  subscript  h  denotes  values  of  the  variable  in  cpiestion  at  the  point  r  =  vi  —  6.  In 
obtaining  Eq.  (7),  we  have  assumed  that  Cr  and.  since  .4  changes  only  by  a  small 
amount  in  any  calculation  of  interest  to  us,  we  have  also  taken  dAldt  and  OA/dz  to  be  zero. 
We  will  henceforth  assume  that  u'j,  =  0  in  accordance  with  previous  assumptions  (I’owell 
and  Zielinski  1959}  that  the  ablated  imoducts  enter  with  no  velocity  and  neglect  the  last 
term  on  the  left-hand  side  of  Eq.  (G). 

In  order  to  obtain  the  rate  of  ablation,  it  is  necessary  to  observe  that  Eqs.  (1)  and  (3), 
which  express  conservation  of  mass  and  energy,  hold  not  only  in  the  plasma  but  iii  the 
solid  polyethylene  as  well  provided  we  take  i;  =  0.  .As  a  result,  the  equations  imply  certain 
boundary  conditions  that  must  hold  at  the  interface  between  the  plasma  and  the  solid. 
These  conditions  can  be  obtained  by  integrating  each  of  the  equations  across  the  boundary 
of  thickness  8  and  then  taking  the  limit  as  6  — ♦  0.  We  finally  observe  that  the  integral  of 
each  partial  time  derivative  ^  produces,  except  for  terms  which  vanish  with  <5.  the  result 
-fj,  <  F  >.  Here  /{,  denotes  the  velocity  with  which  the  boundary  moves  and  the  brackets  (in 
this  section  only)  denote  the  change  in  the  quantity  in  question  as  the  boundary  is  crossed. 
Carrying  out  the  two  integrations  then  and  solving  the  results  for  ut  and  fj,,  we  obtain 


Ub 


Pb  +  PbPb 


(8) 


and 


rb 


Pb-^q 

p,{Pb  +  PbPb) 


(9) 


In  those  two  equations  p,  denotes  the  density  of  solid  polyethylene,  and  A7  =  qrb  —  qra 
denotes  the  difference  in  the  flux  emitted  by  the  gas  at  the  boundary  and  that  transmitted 
all  the  way  into  the  solid  or,  in  other  words,  the  flux  ab.sorbed  by  the  boundary,  In  deriving 
these  equations  wc  have  assumed  that  p,  >>  Pb,  E,  <<  Ei,,  and  ul  «  Eb.  Equation  (8) 
can  evidently  be  used  to  obtain  the  ablation  rate,  whereas  Eq.  (9)  yields  the  rate  at  which 
the  solid  surface  recedes. 


Equations  (8)  and  (9)  can  now  be  used  in  Eqs.  (5)-(7).  We  find  after  expanding  out  the 
derivatives  and  simplifying  the  equations  considerably  that  the  results  can  be  written 


dp  dw 

'dt'^'^'dl 


(10) 
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and 


c)P 

dz 


dll' 


P 


dt 


dw 

+  /’K'—  +  (/•/)„  =  - 


DE  OE  Jw  2  % 

P-^  +  pu’-;^  +  P-;^  +  pa{E  -  uP  2)  - - qr,  -  -7^. 

Ol  dz  uz  a  ri,  uz 


with  the  ablation  rate  pa  given  by 


(12) 


2pb^q 

rb{Eb  +  pbEb)' 


(13) 


The  overbars,  denoting  average  values,  have  for  simplicity  been  omitted  in  these  equations 
and  will  be  omitted  in  the  remainder  of  the  report.  It  is  to  be  understood  that  averaged 
values  are  implied  unless  stated  otherwise.  Equations  similar  to  those  above  have  been 
derived,  for  example,  by  Tidman,  Thio,  CJoldstein,  and  Spicer  (198G). 


The  ablation-dependent  terms  in  the  above  equations  can  be  interpreted  as  follows.  The 
term  pa  in  Eq.  (10)  represents  the  rate  at  which  the  plasma  density  increases  because  of 
ablation.  The  term  wpa  in  Eq.  (11)  represents  a  drag  force  on  the  gas  because  the  ablated 
products  enter  with  no  velocity.  Finally,  in  Eq.  (12),  the  term  paE  accounts  for  energy 
expended  by  the  plasma  to  bring  about  the  ablation  process,  whereas  the  term  Paiv'^/2 
corresponds  to  a  dissipative  effect  resulting  again  from  the  ablation  products  entering  the 
mainstream  of  the  arc  with  no  velocity.  The  first  of  the  last  two  terms  results  in  a  decrease 
in  internal  energy,  while  the  second  results  in  an  increase. 


3.2  Equations  of  State  and  Supplementary  Relations.  As  is  well  known,  the 
con.servation  equations  presented  in  the  previous  section  are  not  sufTicient  to  determine 
uniquely  the  various  fluid  parameters.  Additional  relations  required  include  expressions 
for  the  heat  flux  q  and  conductivity  cr,  a,s  well  as  equations  of  state  w'hich  relate  P  and  E  to 
the  density  p  and  temperature  T.  .Much  of  these  data  have  been  tabulated  for  polyethylene 
in,  for  example,  the  SESAME  tables  (Holian  1984).  Nevertheless,  we  will  formulate  the  ap¬ 
propriate  relations  in  their  entirety  and  carry  out  the  calculations  within  the  computer  code. 
The  puri)ose  for  undertaking  the  calculations  rather  than  relying  on  tabulated  or  curve- 
fitted  values  is  our  interest  in  developing  a  mechanism  whereby  plasma  constituents  may  be 
changed  easily,  even  for  mixtures  whose  properties  have  not  been  previously  documented. 

In  all  the  calculations  to  be  discussed,  we  will  neglect  both  the  longitudinal  heat  flux  q, 
as  well  as  the  flux  conducted  into  the  solid  polyethylene,  qr,.  The  first  effticl  was  deemed 
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negligible  because  the  length  of  our  capillary  lubes  is  generally  large  compared  to  their  radius. 
Moreover,  earlier  calculations  which  we  have  performed  including  this  effect  suggested  onl> 
minor  differences  with  cases  in  which  it  was  neglected.  The  second  effect  is  assumed  to  be 
fairly  insignificant  because  of  the  low  thermal  conductivity  of  polyethylene.  The  energy  flu.x 
radiated  out  the  sides  of  tlie  cairillary  will  be  assumed  to  obey  the  standard  blackbodt- 
radiation  result.  Con.sequcntly,  we  will  neglect  the  last  two  terms  on  the  right-liand  side  of 
Eq.  (12),  and  take  ^<7  in  Eqs.  (8)  and  (9)  to  be  given  by 

=  qrb  =  ersTt,  (M) 


where  as  is  Stefan’s  constant. 

To  determine  the  pressure  equation  of  state,  we  assume  the  law  of  partial  pressures  and 
write 


P  =  ric(l  +  iic  +  2x2c)I^-T  +  n//(l  +  TiH)kT,  (15) 

where  k  is  Boltzmann’s  constant,  where  nc  and  n//  denote  the  number  densities  of  heavy 
particles  for  carbon  and  hydrogen,  respectively,  and  where  is  used  to  denote  the  ratio 
of  "jth”  ions  to  heavy  particles  for  species  a.  We  assume  that  the  carbon  atoms  are  ionized 
at  most  doubly,  an  assumption  which  will  prove  reasonable  for  the  temperature  and  density 
ranges  of  interest.  It  is  sometimes  convenient  to  express  Eq.  (15)  in  terms  of  the  mass  density 
and  to  use  the  rather  obvious  relation 


nc  =  - - ,  (Ifi) 

mac  +  ramoii 

where  niuc  and  mo//  arc  the  masses  of  the  carbon  and  hydrogen  atoms,  re.spectively,  and  ro 
is  the  ratio  of  hydrogen  to  carbon  heavy  particles.  Eor  the  assumed  homogeneous  mixture 
of  polyethylene,  ro  —  2. 

The  ion  concentrations  which  appear  in  Eq.  (15)  can  be  determined  by  solving  the  Saha 
equations  (sec,  c.  g.,  Zcldovich  and  Raizer  19G6).  In  rather  simplified  notation  these  govern¬ 
ing  equations  caii  be  written 


■Tjo  —  A  20  3"  la  — 


h\oh\a _ 

1  d  A  la  -f  A  lo  A  2a 


The  functions  Kja  are  given  by 


(17) 
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2  Z,:  ■2-in,kT..,.,  .  , 

—  y  (  T~,  )  ~  ^  J  \- 

?ie  Ir 


[IS) 


where  Tn,.  is  the  electron  mass,  h  is  Planck's  con.-'taiit.  denotes  thejtli  ionization  potential 
for  specie  q,  and  A/j^  denotes  the  reduction  in  the  ionization  potential  because  of  the 
nonideal  effects  discussed  earlier.  .As  pointed  out,  we  use  the  model  proposed  by  Ebeling 
and  Sandig  (1973)  to  compute  ihe  reduced  potentials.  The  reduction  can  then  be  written 


-lTto(Az,-f  A/S)' 

In  this  expression  Cq  is  the  permittivity  of  free  space,  c  is  the  electronic  charge,  and  .\  is  the 
deDroglie  wavelength  given  by 


h 

{2i:m,kTyl^‘ 


(20) 


The  parameter  Xd  is  the  standard  Debye  length,  taken  here  to  account  foi  both  electron  and 
positive-ion  shielding,  and  can  be  written 


Xo  = 


iokT 


11/2 


[nee2(l  +  Z) 


where  n,  is  the  electron  d<;nsity  given  by 


(21) 


ne  =  nui  +  n,c  +  2n2c, 

and  Z  represents  the  effective  charge  on  an  ion,  namely. 


(22) 


ngXic  +  incx-ic  +  nuxui 
TicXic  +  2ticX2C  +  nffXiH 


(23) 


The  electronic  partition  functions  Zja  which  appear  in  Eq.  (IS)  can  be  written  (Zeldovich 
and  Raizer  1966) 


^:c,  =  Yl(k3u,exp{~UjoilkT),  (2‘1) 

t 

where  t/jo,  is  the  energy  level  of  the  ith  electronic  state  of  the  jth  ion,  and  is  the 
appropriate  degeneracy  factor  for  this  level.  These  energy-level  data  have  been  compiled 
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Table  1.  Degeneracy  Factors  and  Electronic  Energy  Levels  for  Hydrogen,  Carbon,  and 
Their  Ions 


Carbon 

Hydrogen 

■ 

/,c=11.26  ei- 
;  =  0 

l2c=  24.38  cv 

J  =  1 

/3C=47.SC  cv 

7  =  2 

A// =  13. 60  tf 
;  =  0 

ID 

9oCt 

Uoc,{ev) 

g\c\ 

f'ic.(ei') 

g2Ci 

C2c.(ei’) 

SOhf  1 

Uqh,{cv) 

1 

1 

0.0 

2 

0.0 

1 

0.0 

2 

0.0 

2 

3 

0.0020 

a 

0.0079 

9 

6.50 

3 

5 

0.0054 

12 

5.35 

4 

5 

1.27 

m 

ID 

1 

2.69 

B 

D 

5 

4.19 

1 

9 

7.50 

8 

3 

7.70 

9 

15 

7.96 

10 

3 

8,56 

11 

15 

8.66 

12 

3 

8.79 

13 

9 

8.87 

elsewhere  (Moore  1952).  The  summation  in  Eq.  (24)  should  actually  be  extended  over  all 
levels  having  energy  lower  than  the  reduced  ionization  potential.  For  simplicity,  however, 
we  have  extended  the  sum  only  over  levels  whose  energy  was  less  than  about  9  ev.  Our 
calculations  will  indicate  temperatures  within  the  arc  of  at  most  about  3  ev,  so  the  that  the 
exponential  factors  in  the  omitted  terms  are  at  worst  about  5%  of  those  in  the  leading  terms. 
Approximations  of  this  type  are  fairly  common  and  have  been  discussed  elsewhere  (Cambel 
1963).  The  actual  levels  and  degeneracy  factors  used  in  the  calculations  have  been  indicated 
in  Table  1.  It  should  be  noted  that  in  some  Ctises  several  levels  of  comparable  energy  have 
been  grouped  together. 

In  calculating  the  specific  internal  energy  of  the  plasma  there  is  a  certain  arbitrariness 
in  what  one  defines  as  the  zero  level.  In  calculations  which  include  ablation,  however,  it 
is  necessary  to  account  for  the  energy  required  to  vaporize,  dissociate,  and  ionize  the  solid 
polyethylene  ,  ci-s  well  as  to  heat  the  resulting  products  up  to  the  temperature  of  the  plasma 
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core.  One  should  also  account  for  the  <dectroiiic  excitation  energy  in  both  the  neutrals  and 
the  ions,  lliough  this  contribution  is  ficciuently  ncgligil.de.  If  we  take  the  zero  energy  le\el  to 
correspond  to  an  unionized  gas  at  zero  degrees,  the  ajtjrrojrriate  expression  cati  be  written 


13  ,  3  ,  , 

E  =  -[- k'Tnc{  \  +  J"ic  -f-  2  )  4  -k'l  u ]i(  \  -r  xm)  -r  {I ic  —  A/ic  lurUTir + 

P  2  2 


(/i//  —  +  (Ac  +  he  -  ■-^hc  —  ^hc)nc^2C  +  pEv  +  pEp  + 


nc(l  —  trie  -  tr2c)H’oc  +  nctric^bie  +  nc-r 2c''  2C  +  ri//(l  —  Xi;/)U  o//]- 

The  quantities  E^.  atid  Ep  are  the  specific  vaporization  and  dissociation  energies  of  polyethy¬ 
lene,  taken  to  be  given  by  2  hj/g  and  SO  kj/g.  respectively;  the  quantities  are  the 
electronic  excitation  energies  and  can  be  obtained  from  the  various  partition  functions  via 
the  relation 


U'  = 

^  or  ■ 


(26) 


This  representation  of  the  internal  energy  neglects  the  difference  in  the  specific  heats  of  the 
solid  and  gas,  but  the  difference  is  negligible  at  the  high  temperatures  under  consideration. 


Finally,  we  must  obtain  an  expression  for  the  electrical  conductivity  a.  We  account  for 
cdcctron  collisions  with  both  ions  and  neutrals  and  have  from  standard  relations 


—  I  ,  V“  '  / 

where  i/er.  and  Vd  represent  the  collision  frequencies  of  electrons  with  neutrals  and  ions, 
respectively.  If  we  assume  that  the  frequencies  representing  electron-neutral  collisions  can 
be  added,  the  appropriate  overall  frequency  can  be  written  in  terms  the  scattering  cross 
sections  of  neutral  hydrogen  and  carbon  as 


=  l’A(l  -  ^Ic  -  ^2McAc  4-  (1  --  Xi//)nH.4/y], 
where  ly  denotes  the  mean  electron  velocity,  namely. 


V 


(28) 


(29) 


n 


and  where  .-Ic'  and  Aj;  denote  c  ross  seelioiis  of  carbon  and  hydrogi'n.  re.'P<'rti\'ely.  In  all 
calculations,  we  have  used  the  xalues  i  l.oelr  and  Kaplan  Ac  =  30t(7'  and  .1//  =  IT-c* 

wlierc  Go  is  the  Bolir  radius,  approximately  equal  to  5.29  x  10“"  m. 

lor  the  collision  frequency  with  ions,  we  ha\e  used  tiie  Spitzer  formula  l)ut  modific'd  in 
the  niaimer  discussed  carlicT  in  ordiu  to  account  lor  nonideal  effects.  In  particular,  the  well- 
ktiown  "log  .\"  term  in  the  usual  Sjtitzer  expression  is  replaced  by  (Zollweg  and  Liehermann 
1987) 


log  \  ^  log{\  +  (30) 

\vlierc  \m  is  given  by 

r2^c„/,-7' 

-  - 

and  n+  denotes  the  density  of  positive  ions  in  the  plasma.  In  terms  of  other  quantities 
calculated,  we  have 


foA'/ 


rt.c* 


■2/3' 


1  /2 


(31) 


n+  =  nc{Tic  +  X2c)  +  (32) 

As  a  result  of  these  substitutions,  the  Spitzer  collision  frequency  becomes 

3SZn,c'^log{  \  +  1.-1A?^J‘/^ 

I'e.  =  - ™ -  (33) 

where  is  a  fairly  weak  function  of  Z  approximately  equal  to  0.58  at  Z  =  1  and  0.68  at 
Z  =  2.  At  points  between  the  extremes,  the  parameter  can  be  approximated  by  a  linear 
function  with  little  error.  It  is  possible  to  show  that.  Eq.  (33)  reduces  to  the  usual  Spitzer 
result  for  sufficiently  low  densities  and  sufTiciently  high  temperatures. 

3.3  Solution  of  Governing  Equations.  The  formalism  developed  in  the  previous 
two  sub.sections  is  sufli'  ient  to  determine  the  basic  flow  properties  of  the  plasma  in  the 
capillary  arc.  In  particular,  if  we  regard  p,  uq  E,  pa,  T,  and  P  as  the  principal  variables,  then 
Eqs.  (10),  (11),  (12),  (13),  (25).  and  (15)  are  sufficient  to  determine  these  unknowns.  The 
differential  equations  must  be  solved  subject  to  some  set  of  initial  and  boundary  conditions. 
One  obvious  condition  is  that 
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u{:  =  Oj  =  0. 


(31) 


since  the  breech  end  of  the  capillary  is  closed.  There  is  an  additional  coiulilion,  which  is  usu 
ally  taken  to  be  the  flow  rate  or  perhaps  the  pressure,  that  inu.st  be  apidied  at  the  capillary 
exit.  The  additional  condition  arises  because  the  gas  in  the  capillaiy  is  in  general  cou{)led  to 
conditions  outside  the  cajjillary,  and  the  gas  beyond  the  exit  point  i.s  not  considered  in  the 
calculation.  Various  arguments  have  been  made  to  suggest  that  the  flow  is  choked  at  exit 
provided  the  plasma  is  exhausted  into  gas  contained  in  a  sufficiently  large  volume  and  we  will 
generally  apply  some  boundary  condition  that  approximates  this  situation.  If  longitudinal 
heat  conduction  is  accounted  for.  additional  boundary  conditions  must  be  specified  concern¬ 
ing  either  the  temperature  or  heat  flux  at  both  the  breech  and  exit  ends,  but  generally  heat 
transport  in  tlie  axial  direction  has  been  neglected  in  the  calculations.  The  initial  conditions 
on  the  various  variables  will  be  discussed  in  Sec.  G  when  specific  calculations  are  undertaken. 
Finally,  it  necessary  to  assume  values  for  the  flow  variables  just  inside  the  capillary-wall 
boundary,  i.  e.,  Tb,  pb^  Ph,  and  Eb  in  Eqs.  (13)  and  (M).  Unless  otherwise  stated  we  will 
assume  that  these  values  are  approximately  given  by  tneir  values  in  the  mainstream  of  the 
arc.  Such  an  assumption  is  consistent  with  the  one-dimensional  nature  of  the  model. 

In  order  to  solve  the  differential  equations  represented  by  Eqs.  (10)-(12),  we  represented 
the  derivatives  in  the  usual  finite-difference  approximation.  The  resulting  nonlinear  algebraic 
equations  were  solved  by  a  completely  implicit,  modified  Newton-Raphson  method  (Pizer 
1975).  Various  values  of  the  grid  spacing  and  time  step  were  employed,  but  generally  we 
found  that  Ar  =  0.02f  and  A<  =  5  x  10“®s  provided  sufficient  accuracy.  An  artificial 
viscosity  wcis  included  in  the  time-dependent  calculations  to  extend  shock  fronts  over  several 
grid  spacings  if  negative  velocity  gradients  appear  in  the  calculations.  The  specific  expression 
used  was  that  of  von  Neumann  and  Richtmyer  (1950),  viz., 


and 


Q  =  -p 


.  div 


0 


(35) 


Q  =  0,  for  >  0. 

Oz 

In  this  expression,  b  is  an  arbitrary  number  of  order  unity  and  we  have  used  the  value  6  =  1.5. 
The  artificial  viscosity  is  included  in  the  numerical  calculations  by  simply  replacing  P  by 
P  —  Q  in  Eqs.  (11)  and  (12). 
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4.  THE  STEADY-STATE,  ISOTHERMAL  MODEL 


It  will  become  e\  idcnl  when  we  umieitake  some  of  the  mamerical  solvitioiis  of  tlu'  general 
lime-depeiuh’iit  equations,  lliat  th-e  t en.ipei at nre  remain.'  burly  nearly  coirstanl  aeros.s  the 
length  of  the  arc.  h m liieritiori'.  in  most  case.s.  the  current  changes  on  a.  timescale  that  i.s 
slow  compared  to  the  rate  at  which  \aiiuus  hydiods  iiamic  jiiocesses  occur  within  the  plasma. 
I  nJer  such  conditions  it  seems  reasonable  to  ajrprf  imate  the  solution  of  the  cciuations 
in  Sec.  3  with  an  isothermal,  st('ady-state  .solution.  Careful  examination,  of  the  equations 
reveals  that  a  consistent  isothermal  solution  can  be  found  provided  the  following  additional 
assumptions  and  ajrproximat ions  are  made;  (i)  in  addition  to  the  temperature,  the  ion 

icentrations  Tie-  J'2Ci  ^itd  J'l.if.  as  well  as  the  conductivity  a.  are  also  position  indejiendetit ; 
.  d  (ii)  the  kinetic  cncrg\'  of  the  |)iasnta  can  be  neglected  relative  to  its  internal  eiiergv. 
1  hesc  assumptions  will  be  shown  to  be  reasonable  when  tin'  more  general  calculations  are 
undertaken.  The  first  is  acceptable  because  both  the  ionization  state  of  the  plasma  as  well  as 
its  conductivity  arc  sensitive  functions  of  temperature,  but  onlv  weak  functions  of  j)rcssure. 
The  second  assumj)tion  can  be  shown  to  hold  reasonably  well  at  the  high  temperatures  where 
there  exist  many  contributions  to  the  internal  energy.  It  can  furthermore  be  shown  from  the 
governing  equations  that,  when  these  two  assumiJtions  hold,  both  pa  and  E  are  also  constant 
in  s|;)ace. 

If  we  now  make  these  approximations,  write  t!ic  governiiig  differential  eejuations  in  con¬ 
servation  form,  and  integrate  from  zero  to  c.  we  obtain  the  purely  algebraic  results 


pw  =  Pa-  (30) 

pw'  +  P  -  (37) 


and 


pivE  +  Pw  —  J^z/a.  (38) 

In  obtaining  these  equations,  we  have  used  the  boundary  condition  on  ti.'  representc.'d  by 
Eq.  (31)  and  denoted  by  /q  the  pressure  at  the  bre(  ch  end  of  the  capillary.  The  determination 
of  Pq  must  come  from  the  additional  boundary  cc  idition  alluded  to  in  Sec.  3.  3  referring  to, 
say,  the  velocity  of  the  gas  at  the  capillary  exit. 
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If  uc  now  .siiljst il  iilc  I'.q.s,  ('If)),  (l.I).  and  (11)  into  (d^),  wc  caii  sohc  formally  for  tlie 
lciii|ici  at  lire  y  lo  inodiu  o 


.J^r, 

‘2rr/7< 


1/1 


(.•]!)) 


'Jilt'  iiliuvi'  equation  oxprchsns  that  under  the  stated  assumptions  the  energy  dissipated  via 
Joule  healing  is  equivalent  lo  that  radiated  out  the  ea|nllary  walls,  with  the  kinetic  eiierg\' 
Vicing  negligihle. 


U'e  now  ohdei  ve  that  the  jiressure  equation  of  stale  c.ni,  according  to  Kqs.  (15)  and  (16), 
li'’  written 


^'  =  <^',0/'.  (10) 

where  f ',u  repreiiciils  the  ion  aeoustie  speed  in  the  plastna  under  the  assumed  isothermal 
(oinhtioas  and  is  given  explidlly  hv 


a 


1(1 


( 1  -f  J'lr  +  2xjc)kT  4-2(1  +  x\it)kT 
»UoC'  +  ‘•i'no// 


i/J 


('ll) 


It  may  he  ohserved  that  for  the  approNimations  discussed,  (■,{)  is  also  indepcndcMit  of  ju  sition 
ill  the  cnjdllnry,  Huhsliluting  this  relhiionship  for  P  into  Kq,  (37)  and  solving  (3G)  and  (37) 
for  ic  uiid  /■  pK/dilcc  the  explicit  lelat iutishijis 


and 


(12) 


I, (^3) 

M)  ^  j 

Siliiihii  ex|i|e^Nioiis  Im  \r  and  /<  have  heen  ohtahied,  for  example,  hy  Tidtiian,  Thio,  Goldsltdn, 
nnd  Spiiei  (l!)86).  It  is  evident  fiuin  the  resulis  that  n  steady  solution  exists  otily  fur 


J'u  il  d^Pnf'jO’ 


(Id) 


W  ill'll  the  e(|U/ilily  lioldh,  the  vi'loi  jly  at  exit  is  a  maxiimun  and  given  hy  the  isothermal 
holJnd  Hperd  (  ',t). 


Even  with  the  preceding  simplifications,  a  conipioicly  analytic  >oiulioii  tu  the  governing 
equations  is  still  not  possible  since  the  conductivity  and  ion  concentrations  are  very  compli¬ 
cated  functions  of  the  temperature  T.  Since  these  quantities  are  spatially  constant,  however, 
the  solution  to  the  approi)riatG  equations  can  he  found  very  easily  by  iteration  and  the  it¬ 
eration  need  noi  be  carried  out  at  every  grid  point  in  the  domain.  More  si-)ccifically,  vve 
solve  Eqs,  (17),  (27).  and  (39)  for  Jic,  Jjc,  ^in-  T  via  iterative  techniques;  position- 

dependent  variables  such  as  p  and  te,  for  example,  then  follow  directly  from  Eqs.  (42)  and 
(43)  without  any  additional  numerical  solutions.  In  the  expressions  for  a  and  the  ionization 
factors,  it  is  necessary  of  course  to  express  position-dependent  quantities  such  the  electron 
density  n,  in  Eqs.  (18)  and  (21)  by  the  appropriate  value  averaged  over  the  length  of  the 
arc.  This  replacement  is  necessary  in  order  to  accord  with  the  assumption  that  the  ioniza¬ 
tion  factors  and  conductivity  do  not  depend  on  po.sition.  Further  study  of  the  equations 
reveals  that  all  position-dependent  quantities  in  the  relevant  cquation.s  can  be  expressed  as  a 
function  of  the  density  p  and  the  ionization  factors  themselves.  For  example,  one  can  easily 
show  that  Tic  can  be  expressed  as 


^ic  +  2t2c  +  2ti// 

Tlf  —  - - - p. 

rrxac  +  2n7o;/ 


(45) 


An  expression  for  <  p  >,*  in  turn,  follows  by  integrating  Eq.  (42)  fronj  zero  to  (  and  produces 
the  result 


where  S  =  Pal[2if)aC,o). 

5.  EXPERIMENT 

In  thi.s  section  we  describe  the  elements  comprising  the  capillary  assembly  and  the  specific 
measurements  made  on  the  plasma  discharge.  These  measurements  include  the  plasma 
discharge  current,  the  voltage  drop  across  the  capillary  discharge,  and  the  pressure  in  the 
capillary,  all  as  functions  of  time  during  each  of  eight  separate  shots.  Particular  emplicisis  in 
the  di3Cu,ssion  will  be  given  to  the  unique  implementation  of  the  pressure  transducer. 

‘The  br»cj<eti  hcie  and  in  the  remainder  of  the  report  denote  an  averaye  over  the  capillary  length. 
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5.1  Experimental  Set-up.  The  capillary  (liscliargo  ci.a.nriicr  coi.si''t('i,!  of  a  12. 7-;;;?;,' 
outer  diameter.  M2.S-r7?77?  long.  i>olyethyietie  rod.  jires^-fil ted  into  27-7;;?;?  outer  diain'McT 
stainless-steel  lube.  .A  0.35  -  777  77!  diameter,  9o-niru  long,  maraging  steel  rod.  coisted  with 
a  room-temperature  '  ulcanizing  adhesive  tRT\')  and  jjiess-fit  ted  into  the  breech  of  tlie 
capillary,  served  as  the  anode  electrode.  The  end  of  the  electrode  rod  eNtending  from  the 
caj'illary  was  threaded  to  the  positive  bus  from  tlie  power  suppl>'.  The  ciitliode  end  of 
the  capillary  was  joined  to  a  0.25-1-771  long  brass  barrel  with  a  galv;inizcd  pijie  coupling. 
The  joint  was  sealed  by  wrapping  teflon  tape  on  the  pipe  threads  and  compressing  an  0- 
rlng  in  the  interface.  The  capillary  and  the  coii[)ling  were  conpnessed  between  the  two 
C.05  777  X  0.1  m  X  0.26  m  halves  of  an  aluminum  housing.  The  housing  was  bolted  to  the 
ground  plate  on  the  power  supply  and  served  as  the  high-current  return.  The  initial  diameter 
of  the  capillary  was  G.35  vim.  .After  each  firing  the  diameter  of  the  polyethylene  tube  was 
increased  with  the  next  larger  reamer  size.  This  procedure  ])rovided  for  ;i  new  capillary 
diameter  for  each  shot,  but  saved  valuable  time  in  remachining  new  capillary  parts.  The 
length  of  the  ablating  discharge  channel  was  kept  constant  at  60.9  mm. 

The  pulsed  power  supply  used  to  perform  the  experiments  consisted  of  four,  \00-^J 
modules  (Zielinski  1989a,  l9S9b).  Each  2000  fiF  module  used  a  nominal  16  ft  II  pulse-shaping 
inductor  and,  for  near-zero  loads,  provided  for  a  current  pulse  with  a  300  /rs  rise  time  and  a 
transfer  admittance  of  11  kA(kV.  Each  module  had  the  flexibility  to  be  charged  to  different 
initial  voltages  as  well  as  to  be  triggered  independently  in  time.  Ignitrons  (NL-2SSS/\)  were 
used  as  closing  switches  for  each  module  and  diodes  (12-R.A204420’s)  were  included  across 
each  module  to  prevent  voltage  reversal  on  the  capacitor  bank.  .All  module  inductors  were 
connected  to  the  positive  bus  which  served  to  deliver  power  to  the  load.  For  these  shots  the 
power  supply  was  configured  such  that  a  sinusoidal-type  current  pulse  was  j^roduced  having 
a  peak  current  sufficient  to  yield  plcisma  core  temperatures  greater  than  2  cv.  The  jiulse 
width  so  greatly  exceeded  the  duration  of  the  time-dependent  plasma  behavior,  that  the 
latter  Wcis  expected  to  be  negligible.  Typically,  the  initial  capacitor  voltages  were  such  that 
the  first  module,  which  initiated  the  plasma,  weis  charged  to  1  kV ,  and  Modules  2  -  4  were 
each  charged  to  a  progressively  higher  value  of  voltage  (2  -  5  kV).  .A  100  /rs  time  delay  was 
used  between  module  triggers.  A  simple  illustration  of  the  circuit  is  shown  in  Fig.  2. 

A  3S-gage  aluminum  wire  (~  O.lM-mm  diameter),  placed  along  the  axial  length  of  the 
capillary,  was  used  to  initiate  the  plasma  discharge.  The  breech  end  of  the  wire  was  secured 
between  the  anode  rod  and  the  capillary  wail,  and  the  muzzle  end  between  the  capillary  and 
the  breiss  barrel.  The  D.C.  resistance  was  measured  between  the  breech  and  the  aluminum 
housing,  and  was  verified  to  be  less  than  1  fl  before  each  shot  to  ensure  a  reliable  arc 
initiation.  An  insulated  block-and-screw  arrangement  was  secured  behind  the  anode  conncc- 
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tion  to  prevent  loss  of  capillar\'  pressure  and  voltage  breakiiow  ii.  As  the  scicn'.  i  igliicr,e<.l. 
an  axial  force  was  applied  to  help  seal  the  anode  in  the  caiiillary.  A  pieinre  of  the  ;iSse;nl)ied 
capillary,  mounted  in  the  aluminum  housing  and  connc'cted  to  the  ])uUed  ])oucr  supply,  is 
shown  in  Fig.  3. 


Charging  IkxkjU 


Figure  2.  Illustration  of  the  Circuit  for  the  400-kJ  Power  Supply.  Only  the  first  module  is 
shown  explicitly. 


5.2  Diagnostics.  .All  of  the  experimental  electrical  data  were  taken  using  three  Nico- 
Ict  4094C  Digital  Oscilloscopes.  Most  of  the  data  were  recorded  at  a  rate  of  1.0  //s/point 
using  12-bit  amplifiers  (1  MegaSample/s,  amplifier  model  number  45G2).  The  pulsed-powcr 
supply  incorporated  a  single-point  giounding  scheme  with  the  signal  cables  laid  in  close 
proximity  to  the  return  ground  braid  of  the  oscilloscope.  This  positioning  helped  to  avoid 
any  induced  loop  currents  in  the  data-signal  lines.  The  ground  braid  from  the  single- point, 
earth  ground  was  attached  to  the  chassis  of  each  oscilloscope  in  the  control  room,  and  the 
A.C.  conduit  ground  for  each  oscilloscope  was  disconnected.  .At  the  oscilloscope,  each  data 
line  and  each  A.C.  power  cord  wais  looped  several  times  through  a  .separate  ferrite  core.  The 
data  were  stored  on  5j-inch  floppy  disks  and  later  transferred  to  the  large  computer  for 
analysis. 

The  pressure  transducer  used  in  these  experiments  was  a  Kistler  gage,  model  607C.  It 
has  a  range  up  to  689.4  M Pa  and  a  sensitivity  constant  cf  -0.147  pC/psi.  A  radial,  flat- 
bottomed  hole  was  drilled  through  the  stainless  steel  tube  into  the  polyethylene  wall.  .A 
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and  the  capilla,ry  wall.  As  the  firings  progressed,  the  capillary  diameter  was  enlarged  and 


consequently  the  wall  thickness  was  reduced.  A  stainless  steel  mounting  block  was  also  added 


to  the  outside  of  the  stainless  tube  to  [rrovide  extra  support  for  the  mounting  threads  of  the 


gage.  The  gage  was  located  17  nnn  from  the  steel  anode.  .-Vn  access  hole  was  drilled  through 
the  aluminum  housing  with  adequate  clearance  for  insulation  around  the  gage,  the  IlG-174 
miniature  coaxial  cable,  and  the  BNC  connector. 


Figure  3.  Capillary  Experiment  Connected  to  the  400-^’J,  Pulsed-Power  Supjrly. 

Before  each  test  the  gage  was  lightly  coated  with  vacuum  grease,  inserted  in  the  radial 
hole,  and  tightened  while  the  capillary  was  in  a  vacuum  chamber.  This  procedure  reduced 
the  amount  of  air  trapped  between  the  gage  head  and  polyethylene  wall.  Typically,  a  vacuum 
of  a  few  Pascals  (-50  ft  Hg)  was  attained  in  the  chamber  via  use  of  a  roughing  pump. 

During  the  course  of  this  work  consideration  was  given  to  the  acquisition  of  the  pressure- 
transducer  signal  voltage.  Mode;!  607C  is  a  high-impedance,  ballistic  transducer  with  a  rise 
time  of  1 .5  gs  and  is  well  suited  for  dynamic  events.  It  produces  an  output  of  electrical  charge 
proportional  to  its  quartz-crystal  displacement.  Consequently,  a  charge  amplifier  is  needed 
to  convert  the  electrical  charge  into  a  voltage  sigiud  itroportional  to  the  applied  pressure. 
In  these  experiments  a  Kistler  Dual  Mode  Amplifier  (iModcl  5004)  was  u.sed  with  a  short¬ 
time-constant  setting  (i.  e.,  maximum  pre.ssurc  pulse  duration  of  hundreds  of  milliseconds). 
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AddiUonall y.  tlic  o\it[)ut  of  tliis  of  amplifier  (±  10  \’  full  scale)  is  curreiU  limited  and 
may  not  preserve  the  high-frequciiev  content  of  the  lransduc<M'  signal,  ])ai  ticalHrly  when  used 
to  drive  long  coaxial  data  lines.  For  this  reason  full-scale  output  was  kept  around  4  1 

The  large.  rai)idly  time- varying  electric  and  magnetic  fields  near  the  e.xperimeiit  seemed 
to  dictate  the  use  of  a  fibcr-oj)tic  data  line  to  ol>tain  consi.stcnl,  noise-free  pressure  signals. 
.\n  amplitude-modulated  (AM)  transceiver/receiver  pair  (Manage  Incorporated,  model  num¬ 
bers  FOL-100,  PC/Rx,  and  FOL-lOO,  PC5/Tx)  was  employed  initially.  However,  the  gain 
of  tlic  AM  pair  was  found  to  be  both  time  and  amplitude  dependent.  Because  of  tlic  un¬ 
certainties  in  calibration  a  simple,  double-shielded  coaxial  data  line  from  the  cxj)crimcnt 
to  the  oscilloscope  was  tried.  Subsequently,  the  double-shielded  technique  was  employed 
throughout  the  remainder  of  the  experiment. 

The  miniature  coaxial  cable  from  the  pressure  transducer  was  roughly  1  m  long,  and 
located  inside  a  .section  of  ground  braid  (second  shield)  which  was  connected  at  the  stainless- 
steel  capillary  tube  (essentially  at  ground  potential).  The  other  end  of  the  ground  braid  was 
attached  to  a  steel  box  which  contained,  but  was  insulated  from,  the  charge  amplifier.  Signal 
noise  could  easily  be  introduced  into  the  data-line  shields  simply  by  allowing  the  metal  box 
or  outer  braid  to  come  into  direct  contact  with  the  concrete  floor.  Therefore,  the  metal 
box  was  insulated  from  the  floor  with  a  sheet  of  wood  and  mylar.  Before  the  signal  line 
from  the  pressure  gage  was  connected  to  the  charge  amplifier,  the  cable  was  wound  several 
times  through  a  large  ferrite  core.  This  procedure  impeded  any  ground  currents  which  might 
be  induced  in  the  cable  shield.  The  chassis  of  the  charge  amplifier  was  isolated  from  A.C. 
conduit  ground  by  disconnecting  the  third  prong  on  the  power  plug  and  l)y  winding  the  cord 
several  times  through  a  separate  ferrite  core.  Because  the  A.C.  ground  was  lifted,  there 
existed  approximately  50  A.C.  between  the  second  outer  braid  and  ground.  Before  the 
miniature  coaxio.l  cable  could  be  connected  to  the  pressure  transducer,  the  outer  braid  and 
experiment  ground  were  temporarily  connected.  In  this  way  any  potentially  harmful  gage 
current  was  discharged  to  ground.  Once  the  signal  line  w-as  connected,  the  temporary  juniper 
could  be  removed  and  a  more  substantial  ground  connection  could  be  made.  The  aluminum 
housing  w'cis  connected  to  the  single-point  earth  ground  located  at  the  experiment.  The 
output  coaxial  cable  from  the  charge  amplifier  weis  double  shielded  and  connected  at  the 
metal  box  and  at  the  o.scilloscopc  irijiut  BNC  connector.  Careful  attention  was  placed  on 
ensuring  that  the  outer  ground  shield  was  well  insulated  along  the  way  to  the  oscilloscope. 
Before  each  firing,  a  known  input  voltage  was  applied  to  the  charge-amplifier  calibration 
input.  This  procedure  was  followed  to  check  the  charge-amplifier  calibration  rus  well  as  to 
verify  the  amplificr-to-oscilloscope  signal  path. 

Additional  electrical  diagnostics  included  inctisurcmcnts  of  the  load  voltage  and  the  time 
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rates  of  change  of  the  plasma  and  individual  inodulc-iiidurtor  currents.  The  voltage  was 
measured  between  the  anode  connection  block  and  the  aluminum  housing  l>y  using  a  Pearson 
current  transformer  and  resistor  siring.  The  time  rale  of  change  of  llie  currents  was  measured 
v.-ith  Rogowski  coils  placed  around  the  i)Ositive  bus  of  the  pulsed-ptjwor  supj)ly  and  each 
individual  inductor  lead.  These  techniciucs,  used  previously  in  railgun  ex])orimcnls,  are 
described  elsewhere  (Zielinski  and  I’owcll  1990). 

The  plasma  discharge  current  ty|)ically  rose  to  peak  in  about  450  ps.  Monitoring  cacli 
of  four,  separate,  rapidly-varying  curr<;nt  puls<;s,  using  Rogowski  coils,  provided  a  valuable 
diagnostic  since  recorded  dijdt  signals  may  contain  information  on  functional  anomalies 
tliat  an  inlegrated  waveform  would  mask.  The  dijdt  signals  were  recorded  directly  from  the 
Rogowski  coils  and  the  integrations  to  obtain  current  performed  nuiiK.'rically  using  a  five- 
point  smoothing  algorithm.  The  Rogowski  coils  used  for  tliesc  experiments  had  rise  limes  of 
less  than  3  gj.  Tlic  error  between  the  sum  of  the  four  measured  inductor  module  currents 
and  the  measured  total  arc  current  for  all  the  shots  was  typically  less  than  h%. 

5.3  Measured  Data.  There  were  eight  firings  performed  in  this  series  of  experiments. 
Current  and  voltage  data  were  obtained  on  all  firings;  five  firings  produced  valid  pressure 
data.  The  pressure  data  obtained  on  the  remaining  firings  were  not  useful.  Solder,  used 
to  fill  unused  fiber-optic  velocity  measurement  holes  in  the  brass  barrel,  blew  out  owing  to 
capillary  pressure. 

Physical  measurements  in  tliese  experiments  included  the  capillary  diameter  before  and 
after  five  of  the  firings  and  limited  video-tape  footage  during  the  firings.  The  changes  in 
dimensions  of  the  capillary  were  used  to  estimate  the  ablated  mass  and  the  recession  rale 
of  the  capillary  wall.  The  video  data  revealed  the  large  amounts  of  light  and  gas  being 
cxijcllcd  from  the  brass  barrel  many  frames  after  the  depletion  of  the  capacitive  stored  energy. 
Although  only  qualitative,  these  data  provided  insight  into  additional  energy-conversion 
mechanisms  which  occur  on  a  time  scale  greater  than  the  width  of  the  current  puksc. 

The  data  sets  from  each  firing  exhibited  similar  features  and  therefore  only  one  firing, 
Shot  8,  will  be  discussed  in  detail.  .Summaries  of  measured  data  from  the  remaining  shots 
an;  presented  in  Tables  2  and  3  and  will  be  given  more  cursory  discussion  later. 

The  arc  current,  i{i),  was  cornjjrised  of  current  supplied  by  the  four  independently  trig¬ 
gered  capacitor-bank  modules.  For  Sliot  8,  the  initial  charge  voltage.s  were  2.0,  3.5,  4.5,  and 
5.0  kV y  on  Banks  1-4,  respectively.  Between  each  trigger  pulse  were  100  p.s  time  delays. 
The  inductor  currents  of  each  module  arc  .shown  in  Pig.  4  as  a  function  of  time.  The  first 
two  current  pulses  rise  to  their  peak  in  roughly  200  ps.  The  rise  time  for  the  current  of 
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the  rcriiiiiiiing  modules  apjjroiiclie.s  the  short-circuit  value  (300  fis),  owing  to  the  lowered 
resistance  of  the  arc  at  the  higher  current.  'I  he  currents  i)icseiited  in  the  (igui<‘  are  lupolar 
even  though  there  arc  diodes  across  each  module  to  prevent  voltage  re\ersal  on  the  capacitor 
bank.  The  interaction  between  the  cunent -dependent  load  voltage  and  the  module  current;; 
causes  the  capacitor  voltage  to  oscillate  about  a  value  above  zero  [rotential.  I'urtheiinore. 
on  this  shot  it  a|)i)eared  that  the  first  module  had  been  recharged  to  a  higher  vohnge  hv 
the  interacting  currents.  Tliis  cfTect  produced  a  larger  j^cak  current  in  the  first  inductor  at 
t  ~  1100  Its.  Unfortunately,  there  were  not  enough  data  channels  to  monitor  a  comiilete  set 
of  electrical  parameters  for  the  j)ower  supply  and  their  interaction  with  the  plasma  load. 

-  Senk  ' 
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Figure  4.  lour  Measured  Module  Inductor  Currents  a.s  a  Function  of  'fime  for  Shot  8. 

The  total  measured  arc.  current  is  shown  in  Fig.  5.  The  current  waveform  has  a  risc-to- 
peak  of  approximately  450  its.  The  inllcction  points  occurring  every  100  /as  are  indicative  of 
the  triggered  niodiilec  supplying  current  to  tlie  load.  'J'hc  final  pulse  width  i.s  about  8.50  jis. 
The  small  modulation  in  the  load  current  for  early  limes  has  been  observi.-d  jrreviously 
during  power-supply  characlcrizutioii  tests  with  fixed-load  impedances.  The  amplitude  of 
tlic  modulations  can  be  controlled  by  adjusting  the  load  impedance.  Exact  calculation  of 
this  effect  is  imjiossible  inaiidy  because  of  the  undeterminable,  distributed  stray  capacitances 
and  inductances  in  the  pulsed- posver  supply.  However,  an  excellent  description  of  such 
phenoinena  and  their  relation  to  measured  transient  data  can  be  found  elsewhere  (Hurdcri 
and  Shear  19G9).  In  summary,  it  has  been  found  when  the  load  impedance  is  large  compared 
to  the  module  impedance,  the  load  current  tends  to  luive  many,  small-aiii|)litude  oscillations 
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superimposed  on  the  load-current  \va\cforni.  Such  an  effect  also  occurs,  as  seen  in  I'ig.  5,  for 
times  less  than  40  fis,  when  the  plasma  arc  is  initiated.  When  a  relati\el_v  small  impedance 
is  used  for  the  load,  as  is  the  case  when  the  current  rises  to  a  large  level  (t  >  200  /ns),  very 
few,  smaller-amjdil  ude  oscillations  apjx'ar. 
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Figure  5.  Measured  Plasma  Arc  Current  as  a  Function  of  Time  for  Shot  8. 

Figure  G  shows  a  trace  of  the  measured  voltage  for  Shot  8.  This  voltage  includes  the 
capillary  resistive  and  inductive  t(;rms.  The  first  18  //s,  during  which  peak  voltage  reaches 
4.5  kV,  corresponds  to  vaporization  of  the  aluminum  wire.  For  times  less  than  100  fis  and 
greater  than  18  /iS,  the  voltage  waveform  reveals  several  sriiail  oscillations.  These  oscillations 
occur  for  the  same  reasons  as  those  discussed  for  the  arc  current  in  Fig.  4  and,  to  a  lesser 
extent,  because  of  the  early,  time-dependent  formation  of  the  plasma.  For  the  remaining 
pulse  duration,  the  voltage  drop  roughly  follows  the  shape  of  the  current  pulse,  with  some  fine 
structure  appearing  in  the  waveform  after  peak  current.  The  capillary  voltage  drop,  V',  was 
obtained  by  subtracting  the  capillary  inductive  term  Lenp^^^  and  the  steel-anode  resistive; 
term,  Ranode^it) ,  from  the  measured  ariode-to  housing  voltage  drop  as  a  function  of  time. 
Impedance-bridge  measurements  of  the  cajjillary  inductance  [Leap)  and  anode  resistance 
{Ranoii)  w'cre  made  at  a  freciuciicy  of  1  kHz,  and  are  101.5  nil  and  1.75  mfl,  respectively. 
The  inductive-voltage  term  contributes  less  than  3%  to  the  total  measured  voltage  and  is 
largest  at  early  times.  The  resistive  contribution  is  slightly  larger,  about  5%,  and  is  greatest 
at  peak  current. 

Figure  7  disi)lay8  the  pressure  trace  for  Shot  8  as  a  function  of  time.  It  should  be  recalled 
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Figure  6.  Mccisurcd  Aiiode-to-llousing  Voltage  Drop  as  a  Function  of  Time  for  Shot  8. 

that  this  pressure,  denoted  by  /p,  wa.s  measured  at  a  point  17  mn?  from  the  anode  end  of 
the  capillary.  The  pressure  waveform  is  similar  in  shai)e  to  the  current  waveform  but  reaches 
its  peak  slightly  later  in  time.  Again,  the  inflection  points  which  occur  at  100  ^ls  intervals 
are  indicative  of  current  supplied  to  the  plasma  load  by  the  three  remaining  modules.  Closer 
inspection  of  the  data  reveals  a  slight  time  delay  between  the  current  and  pressure  inflection 
points.  This  difTercnce  results  from  several  effects.  Foremost,  the  gage  has  a  published 
response  of  1.5  /is  and  the  thin  polyethylene  wall  between  the  gage  head  and  capillary  wall 
will  no  doubt  increase  that  value.  Furthermore,  the  plasma  may  exhibit  some  time-dependent 
behavior,  especially  with  rapidly  varying  dildt  in  the  circuit.  For  this  shot  (Shot  8)  the  time 
delay  between  the  current  and  i)ressure  inflection  points  was  4  /is  for  early  times  but  on 
Shot  4,  which  had  a  much  smaller  peak  current  and  pressure,  a  10  delay  was  observed. 
These  value.s  are  consistent  with  the  effects  mentioned  above. 

Ground  potential,  and  therefore  tlic  reference  for  the  oscilloscope,  can  shift  due  to  small 
currents  flowing  in  the  inductance  and  resistance  of  the  ground  connection.  As  a  consequence, 
pulsed  signals  can  converge  to  nonzero  values  for  long  times.  In  order  to  verify  that  the 
pressure  signal  returned  to  the  j)re-shot  zero  baseline,  we  extended  the  measurement  for  one 
shot  (Shot  4)  long  after  the  pulse-width  of  the  current.  In  this  measurement  the  transient 
recorder  acquired  nearly  8000  /ns  of  data  and  the  pressure  waveform  was  found  to  approach 
to  zero  at  8000  /rs. 
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Figure  7.  Measured  Pressure  as  a  Function  of  Time  for  Sliot  8. 

Perhaps,  however,  the  most  interesting  regime  for  examining  both  the  accuracy  and 
acquisition  of  the  pressure  data  is  the  initiation  phase  of  the  plasma  discharge.  The  fine 
cylindrical  aluminum  wire  used  to  initiate  the  plasma  is  exploded  by  the  passage  of  the 
current  pulse  produced  by  Bank  1.  This  method  of  plasma  initiation  is  similar  to  experiments 
performed  previously  utilizing  cxploding-wire  phenomena.  With  wire  diameters  on  the  order 
of  0.1  mm  and  available  energies  larger  than  5  J/cm,  a  nearly  cylindrical  shock  wave  system  is 
produced  (Dennett  1959).  The  system  is  fairly  complex;  however  strong-shock,  similarity  flow 
thcor}  (Lin  195-1)  has  been  found  to  describe  the  trajectory  of  the  experimental  shcck  wave 
and  can  be  used  to  estimate  the  peak  jjressure  when  the  plasma  is  initiated.  In  the  strong- 
shock  approximation  it  is  assumed  that  the  energy  is  deposited  instantaneously.  Uniform 
atmosphere  surrounding  the  wire  is  assumed  to  consist  of  perfect  gas  with  constant  heat 
capacity.  In  our  experiments  energy  deposition  is  not  instantaneous  and,  in  fact,  extends 
over  a  period  an  order  of  magnitude  longer  in  duration  than  that  observed  in  exploding  wire 
experiments  (Burden  and  Shear  19G9).  Therefore,  similarity  flow  theory  should  yield  an 
upper  limit  for  the  pressure  produced  in  the  capillary  during  arc  initiation  and  is  given  by, 

0.216  i:., 

‘sas]  —  r  2 

The  energy  deposited  in  the  conducting  wire,  is  the  value  taken  when  the  capillary  voltage 
reaches  a  maximum.  At  this  time  the  wire  conductivity  is  at  a  mirurnum.  Thereafter,  the 
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air  becomes  increasingly  ionii^cd  and  rcdalively  more  energy  is  dc])osite(l  ir.  ti,c  colr.nin  of  air 
surrounding  the  remaining  metal  wire. 

The  data  for  Shots  4.  6,  7,  and  S  were  used  to  calculate  the  pressure  after  wire  vapor¬ 
ization.  In  all  tests  the  calculated  pres.sure  overestimated  the  mea.sured  pres.-ure  by  a  factor 
of  about  5  at  low  energies  (less  than  1  J/cjii)  and  about  2  at  the  larger  dej)o^ited  energies 
(~  3  Jlcm).  The  recorded  pressure  pulse  from  wire  vaporization  peaks  5-10  /a-  after  the  wire 
exploded.  This  observation  is  consistent  with  the  development  of  a  cylindrical  .shock  wave. 
However,  the  observed  pulse  duration  is  an  order  of  magnitude  longer  than  strong-shock 
theory  would  suggest.  The  differences  in  both  time  and  ntagnitude  between  the  calculated 
and  measured  initiation  event  are  largely  attributable  to  the  gage  rcsjronse  (see  discussion  in 
the  previous  paragraphs)  and  calibration  (sec  discussion  in  Sec.  5.  4).  Further  differences  can 
be  accounted  for  by  deposition  of  nonconstant  wire  energy  in  a  finite  time.  Most  important, 
is  the  absence  in  the  pressure  signal  of  noise  induced  from  capacitor  bank  trigger  pulses  and 
circuit  transients  during  wire  explosion.  This  final  point  gave  us  confidence  that  tiie  signal 
recorded  weis  an  accurate  representation  of  the  transducer  output. 

A  summary  of  the  experimentally  measured  parameters  for  all  the  shots  is  shown  in 
Table  2.  Included  are  the  current,  tfC,  in  Column  2;  the  resulting  values  for  the  capillary 
voltage,  V,  with  corrections  for  the  inductive  and  resistive  terms  in  Column  3;  the  value  of 
the  capillary  pressure,  /p,  in  Column  4;  and  the  time,  at  which  these  data  are  obtained 
in  Column  5.  The  next  two  columns  show  the  value  of  the  peak  pressure.  Ppk,  and  the  time 
at  which  peak  pressure  occurs,  respectively.  The  last  two  columns  list  the  measured 
values  for  the  initial  capillary  radius,  rj,,,  and  the  final  capillary  radius,  ?•(,;. 

Table  3  summarizes  additional,  relevant  parameters  obtained  from  the  measured  electrical 
data.  It  should  be  noted  that  there  appears  to  be  a  time  correlation  among  the  peaks  of 
the  voltage,  power,  current,  and  pressure.  The  peak  voltage  across  the  capillary  occurs  first, 
followed  by  the  peak  in  the  arc  current  roughly  50-75  ps  later.  Because  of  the  oscillatory 
nature  of  the  discharge,  peak  power  occurs  a  few  microseconds  before  peak  current.  Peak 
pressure  follows  peak  current  by  approximately  100  fis.  Contrary  to  common  assumptions 
that  capillaries  of  this  type  operate  in  the  steady-state  regime,  the  measured  ])lasma  data 
appear  to  exhibit  a  substantial  amount  of  time  dependence.  This  point  will  be  further 
discussed  in  the  following  pages. 

Table  4  presents  a  summary  of  results  calculated  from  various  measured  electrical  and 
physical  data.  The  change  in  capillary  radius,  Ar,  listed  in  Column  2.  was  obtained  from  the 
difference  of  the  final  and  initial  values  for  the  capillary  radii  (see  Table  2).  The  capillary 
mass  loss  was  calculated  from  that  difference  and  the  density  of  polyethylene  (950  kglm^). 
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Table  2.  Suinmai  v  of  Measured  Data 
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I\.{MPa] 

r,i,[MPa) 
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i 

4 

31.9 
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91.6 
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3.365 

3.4.54 

5 

41.9 

1 777.4 

473 

141.6 
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3.543 

6 
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1 16.0 
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3.568 
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8 
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HB 

143.9 
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3.835 

t  Not  measured. 

*  Pressure  nieasuremeiits  not  valid, 
t  Soldered  fiber-optic  holes  had  blown  out. 


Table  3.  Peak  Data 
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8 
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120.0 
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*  Pressure  measurements  not  valid. 
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Table  4.  Summary  of  Energy-Conversion  Data 


£D[kJ) 

7  (%^| 
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DH 
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23.4 1 

5.3 

224. 9 1 
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45.3 

1  2 

27.5t 

32.3t 

7.3 

225.3t 
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68.6 

3 

56. dt 

66. ot 

15.0 

224.7t 

62.2 

82.6 

4 

SS.5 

109.7 

21.4 

194.4 

49.8 

81.1 

178.0 

113.2 

29.8 

263.8 

41.7 

83.1 

203.0 

1.33.5 

34.1 

255.4 

40.1 

72.2 

22S.0 

155.0 

37.4 

241.4 

41.7 

71.9 

8 

254.0 

ISO.O 

47.0 

261.-1 

_ 

34.0 

7.3  6 

t  Extrapolated  from  results  on  Shots  4-8. 

The  total  energy  deposited  in  the  capillary,  Eq  (Column  4),  was  obtaip.ed  by  integrating  the 
product  of  the  measured  capillary  voltage  and  current  over  the  width  of  the  current  pulse. 
The  ablation  coefficient  for  the  polyethylene  capillary  was  taken  to  be  the  ratio  of  deposited 
energy  to  mass  loss  and  is  shown  in  Column  5. 

We  should  finally  point  out  that,  for  a  power  supply  delivering  energy  to  a  load,  the 
load  impedance  plays  a  crucial  role  in  determining  the  amount  of  energy  transferred.  In 
electrothermal  guns,  we  are  concerned  with  maximizing  the  energy  dissipated  in  the  plasma, 
Typically,  an  impedance  can  be  found  which  will  maximize  the  amount  of  energy  tran.sferred 
and  hence  the  trai.sfcr  efficiency  t}.  As  the  load  impedance  excecjds  that  for  inaximurn 
efficiency,  the  amount  of  transferred  energy  steadily  decreases;  for  load  in)j)(;danc(;.s  even 
slightly  less  than  that  at  maximum  efficiency,  the  amount  of  energy  transferred  can  be 
greatly  decreased.  This  effect  is  particularly  pronounced  with  a  plasma  load  sinctr  the  i)la.sina 
resistance  becomes  smaller  for  larger  currents.  Even  in  the  relatively  short  serie.s  of  firings 
discussed  here,  the  effect  of  load  resistance  on  the  transfer  efficiency  can  be  seen,  'J'he 
minimum  of  for  each  shot  was  used  to  calculate  the  minimum  plasma  re.sislanc(,’, 

Hrnin^  listcd  in  Column  6  of  Table  4.  Low-current  shots  i)roducc  a  fairly  high-re.sistance 
load  for  which  the  transfer  efficiency  is  found  to  be  less  than  70%  (Shoti)  1  and  2),  I'bi  a 
load  resistance  below  40  r/Tl,  the  transfer  efficiency  also  begins  to  decrease.  For  thi.s  power 
supply,  maximum  efficiency  for  energy  transfer  is  near  8.5%  and  occurs  fur  a  load  impedance 
in  the  range  of  40-65  mfl. 
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b.'i  PrL’Bwurt-GtiKO  Cnlihratiun.  In  our  rally  attriiipls  to  iiicasuic  tlic  i)laMna  ])rr,s- 
hiM:'.  a  Niiinll,  silit  (Uia  grca.^c  lillcd  ludc  w.is  used  to  cuiiiiiiiiiiicate  ijressutc  in  the  waluinc  of 
lln*  fajiillnry  lo  llie  pressure  iransdiner,  lluwever,  I  Ids  idea  was  quickly  adaiuloned  as  the 
iilidilic  gicasc  failed  t(  insulate  llie  gage  from  llie  liigli.  pulsed  electrical  current,  ,\n  altei- 
iinii'.'c  iiiclliud  Wan  Id  simply  insert  the  gage  under  v.icuiiiii  into  a  shallo'.v  cavity  macldned 
in  the  uulei  capillnty  wall,  llieieby  leaving  a  tliin  section  of  polyethylene  between  the  gage 
head  and  the  iniiet  surface  of  the  capillary  wall.  ']'he  transducer  had  a  flat,  circular  steel 
iliapliragiii  weld''<l  to  the  Iransflm'er  rasing.  When  nmnnied  radially,  the  polyethyh’iie  wall 
llieii  pioji'i  led  a  lioiiiiidhji  III  (  u>ss  sei  tion  !(<  the  gage  liiMid,  (me  which  was  thin  <il  the  center 
of  the  gage  and  thicker  at  the  ends, 

lo  esahiaie  the  liniisdmei  system  as  (oinpletely  as  iiossihle,  it  was  decided  to  jierform 
t  he  giige  I  alilii at  mil  in  si/ a  ‘J  o  ohi idn  pteNsiiren  c(»m|iarahle  lo  the  eleriric-discharge  exper¬ 
iments.  a  fluid  and  piston  ariangement  was  employed.  A  short  section  of  l)arrel  fabricated 
foi  (alihialioii  purposes  was  joined  Irr  the  capillary  using  a  galvanized  |dpc  coupling.  A 
>i'i  olid  Isislh'i  |i|essiiie  liali  dm  ei  (Mt/del  useil  foi  a  refeieiice  iiieasiireinent,  was 

limiinled  ill  the  standaid  haliislii  oi  eniaiioii  ne,ti  1  he  muzzde  end  of  the  ralihration  barrel, 
,A  imlye'hyleiie  fod,  drilled  to  the  depth  of  OU.i)  mm,  w;u  used  lo  form  the  capillary.  A  steel 
|dstoii  wlili  lliiee  0  rings  around  llm  piston  slntfl  was  inoiiiiled  on  four  guide  pins  and  was 
free  lo  slide  Into  the  inu/.ji'.le  eml  of  the  calihralion  hariol.  DifruHion  pump  oil  was  used  as 
llie  I uinpieising  Ihiid.  A  stand  was  coiistiucled  by  which  steed  disks  could  be  dro|)|>e(l  onto 
llie  il(i<d  piston,  'I  lie  Weights  (.'10  kfi  inaxiitiurii),  which  we-e  attached  to  an  nliimimim  plate, 
we)e  guided  by  foin  sleid  wiies,  The  (  hpilhiiy  and  (alihralion  barred  were  asseiiilded  in  the 
Aluiiilimiii  hulking  and  p<jjili(»iied  iipiighl  at  l|i«‘  hotloin  of  the  drop  stand,  do  ensure  lliul 
iniiilmal  all  was  Impped  in  the  (^uinpiessing  fluid  pri(;r  lo  insertion  of  the  piston  and  drop¬ 
ping  the  weight,  the  piston  had  a  small  threadeil  hide  through  the  shaft  by  wdiicli  lra))ped 
all  (ould  be  released,  Uin c  llie  air  was  suflitienlly  removed  from  the  barred  and  the  piston 
leHted  lu  the  propel  depth,  the  through  Inde  was  sealed,  'I'lie  wcighis  Were  released  one  me- 
l«  1  iihose  the  piston,  siihsecjiienlly  ( i/iiipiessed  tin*  lliild  III  the  capillary  and  the  ealihialion 
ballel,  and  llieiehy  geiielaled  the  deslied  piessiiie  pulse,  I'euk  pressures  similar  to  those  in 
llie  eli  ilinal  djsLiiaiges  We|e  ohlnliied,  allliougli  pressure  rise  limes  were  on  the  order  of  a 
fev^  iiiilllsci oiids.  'Ibis  iyp»  of  meastiicmciii  was  lepealeil  foi  each  capillary  diameter  fired, 
and  I'll  saiiiiiis  siz' »  of  diop  weight,  A  typi<  al  set  of  picssuie  traces  ici  slncwii  in  big,  8,  It 
(all  hi  seen  (loiii  I  Ig  8  llial  both  ll  aiisilm  ei  signals  (loss  at  early  limes,  'J'his  Indiavioi  is 
paillally  due  lo  the  piopagatiuii  d(day  time  of  the  pressure  wave  In  the  Ihild  cohillin, 

'I  lie  I  III  lo  of  I  III'  piessiiie  |ei  o|i|ed  wil  h  the  lefi'lem  «•  gage  to  llie  piesMlJK*  re(  onled  b)  ( he 
iMplIlai)  gage  Was  taken  lo  he  the  i  alihialioii  fai  loi  and,  as  liidiiated,  was  deleimilied  for 


Figure  8.  Reference  Pressure  and  Capillary  Pressure  as  a  Function  of  Time. 

cacli  capillary  diameter  tested.  The  tests  which  yielded  the  largest  peak  pressures  were  used 
in  determining  the  calibration  factor.  Similar  calibration  curves  were  obtained  regardless  of 
the  j)eak  pressure  for  a  given  capillary  diameter.  The  calibration  curve.s  can  be  separated  into 
two  distinct  regions:  positive  dl'jdi  (rising  pressure)  and  uogal'we  dP/dt  (falling  pressure). 
Curves  for  the  falling  pressure  data  are  shown  in  Pig.  9.  The  calibration  factor  is  fairly 
constant,  ()articularly  for  j)ressurcs  greater  than  the  yield  strength  of  polyethylene  (31  M Pa). 
Furtliermore,  deviation  from  the  reference  traii.sduccr  reading  is  on  the  order  of  a  few  i)ercent 
for  the  larger  diameters  and  15%  for  the  smaller  capillary  diameters.  The  rising  pressure 
data  are  shown  in  F’ig.  lU.  'I’liese  data  exhibit  the  same  behavior,  except  for  some  oscillations 
occurring  at  early  limes.  The  nature  of  this  behavior  can  be  attributed  to  the  inadequate 
seals  and  Joints  used  in  the  barrel  a.ssembly.  During  the  initial  pressure  rise,  the  fluid  probably 
(illetl  voids  funned  during  ao.semhly  of  the  caj)illary  to  the  barrel. 

Of  particular  interest  is  the  calibration  test  for  the  (i. 731  rnrn  diarueter  capillary.  On  that 
le.Ht  the  droj)  weight  hounced  on  the  jjiston  producing  a  daublc-i)eak<!d  pre.ssurc  pulse  for 
wliicli  the  second  pulse  wa.s  iiiitiat(,‘(l  before  the  first  couhl  fall  to  zero.  It  appears  from  these 
data  (hat  once  the  itiilial  pre.s.surization  of  the  fluid  in  the  capillary  coluimi  was  completed 
(at  80  Mj’u)y  the  calibration  factor  remained  fairly  constant  at  0.830  with  a  deviation  of 
±  3.'1%,  even  for  pressures  wliich  subseriueiilly  fell  below  80  M Pa. 

A  calibration  test  was  also  done  using  the  l.V.il-rnin  diameter  capillary;  as  an  e.xpcdiciicy. 
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Figure  9.  Calibration  Factor  as  a  Function  of  Measured  Capillary  Pressure.  Data  are  for 
falling  pressure. 
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Figure  10.  Calibration  Factor  os  a  Function  of  Measured  Capillary  Pressure.  Data  are  for 
rising  pressure. 
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Table  5.  S  ummar\’  of  Calibration  Constants  for  the  Capillary  Gage 


Capillary  Diameter  (mm) 

Wall  Thickness  o  the  Gage  (mm) 

Calibration  Factor 

6.350 

1.77 

0.863 

6.731 

l.oS 

0.823 

6.908 

1.49 

1.044 

7.137 

1.38 

1.035 

the  capillary  gage  was  installed  under  atmospheric  conditions  to  facilitate  the  tests.  There 
was  no  appreciable  deviation  in  the  calibration  factor  for  pressures  greater  than  75  MPa. 
However,  the  oscillations  appearing  at  low,  rising  pressure  were  somewhat  erratic. 

Data  for  the  largest  capillary  diameter  tested  (7.67  mm)  are  not  shown.  The  peak 
pressure  attained  using  that  diameter  (ISO  MPa)  was  sufficient  to  burst  the  galvanized  pipe 
coupling  that  joined  the  capillary  to  the  calibration  barrel.  However,  the  calibration  factor 
obtained  prior  to  failure  did  not  deviate  significantly  from  that  obtained  using  the  6.908  and 
7.137  mm  capillary  diameters. 

It  should  be  noted  that  the  form  and  duration  of  the  pressure  pulses  are  not  identical  to 
those  generated  in  the  electrical-discharge  experiments.  In  the  electrical  firings,  one  end  of 
the  capillary  was  open  and  the  other  sealed  with  the  steel  anode.  In  the  calibration  experi¬ 
ments,  both  ends  were,  to  a  great  extent,  closed.  Also,  the  calibration  tests  were  performed 
at  room  temperature  whereas  the  electrical  tests  produced  elevated  wall  temperatures.  High 
temperatures  should  cause  the  polyethylene  to  yield  at  a  reduced  pressure  allowing  the  gage 
to  follow  the  plasma  pressure  more  closely  as  a  function  of  time. 

The  capillary  gage  calibration  factors  were  taken  to  be  a  best-fit,  constant-value  for 
pressures  greater  than  the  yield  strength  of  polyethylene  (34  MPa),  .^n  average  of  the 
rising  and  falling  data  was  used,  and  are  listed  in  Table  5.  As  can  be  seen  from  the  table, 
wall  thicknesses  greater  than  1.58  mm  cause  an  overestimate  in  the  actual  pressure  by  about 
15%,  while  thinner  wall  thicknesses  yield  underestimated  pressures  by  only  about  4%. 

6.  CALCULATIONS  AND  COMPARISON  WITH  EXPERIMENT 

We  now  turn  to  undertaking  some  specific  calculations  with  both  the  steady-state  and 
time-dependent  models.  In  all  cases  we  attempt  to  carry  out  calculations  which  allow  us 
to  compare  the  results  with  experimental  data.  Considerable  discussion  will  be  presented 
indicating  the  extent  to  which  the  theory  and  experiment  agree,  as  well  as  attempting  to 
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interpret  physically  the  results. 


6.1  Steadv-State  Calculations.  In  PZl  we  exercised  the  steadv-statc.  isothermal 
model  to  compute  j)lasina  properties  for  five  of  the  shots  discussed  in  the  previous  section. 
The  shots  selected  for  analysis  were  those  for  which  valid  pressure  data  were  obtained,  i.  e.. 
Shots  2,  4,  5,  C,  and  8.  The  calculations  were  carried  out  using  as  input  the  peak  current 
obtained  in  each  of  the  five  shots.  A  comparison  of  the  theoretically  obtained  voltage  drop, 
Vth.  measured  across  the  length  of  the  capillary,  was  also  compared  with  that  obtained 
experimentally  and  given  here  by  lir-  The  voltage  can  be  calculated  from  the  simple  and 
obvious  formula 


Vth  =  til(-rla).  (47) 

We  now  undertake  similar  calculations  using  the  updated  steady-state  model.  The  prin¬ 
cipal  difference  between  these  calculations  and  those  reported  in  PZl  is  that  the  present 
model  includes  the  nonideal  effects  discussed  in  Sec.  2.  Results  of  the  calculations  are  shown 
in  Table  6.  Most  of  the  results  are  self-explanatory  and  need  little  discussion.  The  radius  rj, 
was  obtained  from  Ihe  average  of  ri,,  and  rt/  indicated  in  Table  2.  The  capillary  length,  f, 
Wcis  in  every  case  equal  to  6.09  cm. 

The  quantity  Vo,  shown  in  Column  7,  represents  the  power  dissipated  in  the  capillary 
due  to  the  finite  resistance  and  is  given  by  Vo  =  ipkVth-  This  quantity  is  important  in 
a.sscssing  the  capillary  performance  because,  in  the  steady  state  with  no  heat  loss  to  the 
walls,  this  energy  is  equivalent  to  that  which  flows  past  the  capillary  exit."  Thus,  typically, 
for  time  scales  of  the  order  of  a  few  milliseconds,  one  can  use  the  capillary  to  furnish  energy 
of  the  order  of  several  tens  of  kilojoules.  The  quantity  Z  in  Column  8  represents  the  ratio  of 
the  total  number  of  ions  in  the  plasma  to  the  total  number  of  heavy  particles.  Consequently, 
this  number  gives  some  indication  of  the  extent  of  ionization  and  increases  with  increasing 
current  and  temperature  (Column  G).  The  parameter  it,  in  Column  9  gives  the  recession  rate 
of  the  solid  surface  and  m,  in  Column  10,  the  steady-state  meiss  contained  in  the  capillary, 
Finally,  <  P  >  represents  the  pressure  averaged  over  the  capillary  length,  and  <  7^  > 
represents  the  similarly  averaged  nonideal  parameter  which  measures  roughly  the  ratio  of 
kinetic  to  potential  energy  in  the  plasma,  Various,  nearly  similar,  definitions  of  the  nonideal 
parameter  exist  and  we  have  used  the  relation 

*1/  ihe  energy  passing  the  capillary  exit  per  unit  time  is  calculated  in  the  isothermal,  steady-stale  model,  it  will  be  found 
to  differ  from  the  dissipated  energy  by  a  slight  amount  equal  to  the  kinetic  energy  of  the  gas  products.  The  reason,  of  course, 
is  thst  the  kinetic  energy  was  neglected  relative  to  the  internal  energy  in  the  approximations  leading  to  the  development  of  the 
model.  We  are  grateful  to  Dr.  Jad  H.  Balleh  who  pointed  out  to  us  the  existence  of  and  reason  for  this  discrepancy. 
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Table  6.  Calculations  with  Isolliernial,  Steady-Stale  Model  and  Comparison  with  Experi¬ 
mental  Results 


.  Shot 

# 

h'k 

ikA) 

rb 

(n?77!) 

V'.r 

(volts) 

v;n 

(volts) 

T 

(cr) 

Vd 

(.MW) 

■ 

rb 

(cm/s) 

7?i 

<  P> 

(.MPa) 

<  > 

2 

10.5 

3.18 

1440 

1339 

1.83 

14.1 

0.39 

3.74 

0.68 

18.4 

0.26 

4 

31.9 

3.41 

1598 

192G 

2.61 

61.4 

0.71 

9.64 

1.43 

58.8 

0.27 

5 

41.9 

3.50 

1777 

2113 

2.84 

88.5 

0.76 

12.4 

1.78 

78.3 

0.27 

6 

45.3 

3.62 

1832 

2109 

2.87 

95.5 

0.77 

12.8 

1.88 

78.5 

0.27 

1  8 

58.7 

3.90 

2035 

218S 

3.03 

128.4 

0.81 

15.0 

2.28 

89.0 

0.26 

_  c^jn,  + 

"  i-rtokT 


(4S) 


As  WcLS  pointed  out  previously,  values  of  this  parameter  larger  than  about  0.1  indicate  that 
Coulomb  interactions  and  the  resulting  nonideal  effects  should  not  be  neglected. 

Comparing  the  voltages  Ver  and  Vth  for  the  various  shots,  we  see  that  the  agreement 
is  fairly  good;  maximum  deviations  are  about  20%.  We  should  also  point  out  that  the 
results  are  not  vastly  different  from  those  obtained  via  the  model  in  PZl.  Voltages  are 
slightly  higher  for  the  compulations  liere  since  the  nonideal  effects  lead  to  a  reduction  in 
the  electrical  conductivity.  Experimental  voltages  are  slightly  smaller  than  those  in  PZl 
because  corrections  for  anode  resistance  and  capillary  inductance  have  now  been  included 
as  discussed  in  Sec.  5.  For  none  of  the  quantities  computed  in  the  table,  however,  did  we 
observe  differences  larger  than  about  10%  from  those  calculated  in  PZl.  This  conclusion 
is  perhaps  surprising  since,  at  a  given  temperature  and  pressure,  nonideal  effects  can  have 
a  much  more  significant  impact  on  the  conductivity.  In  the  problem  at  hand,  however,  in 
which  all  the  variables  are  nonlinearly  coupled,  a  reduction  in  the  conductivity  because  of 
nonideal  interactions  leads  to  an  increase  in  resistive  heating.  This  increaise  is  followed  by  a 
rise  in  temperature  which  tends  to  bring  the  conductivity  back  up.  Thus,  the  overall  effect 
is  not  so  significant  as  might  be  expected  a  priori. 


In  PZl  we  also  investigated  a  model  which  was  steady  state,  but  in  which  the  isothermal 
approximation  was  not  assumed.  In  other  words,  the  variables  which  were  assumed  to  be 
position  independent  in  the  previous  section,  were  allowed  to  vary.  In  the  present  work,  we 
have  not  developed  such  a  model  explicitly,  but  rather  have  obtained  steady-state  results  by 
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relaxing  the  time-dt‘j)cndent  calculations  for  a  specific  current  until  nearly  steady  bchaiior 
was  observed.  We  discuss  now  some  of  the  results  of  a  calculation  undertaken  for  Sliot  S  (see 
Table  6).  This  calculation  was  performed  primarily  to  examine  the  validity  of  the  isothermal 
approximation. 

The  initial  conditions  for  the  calculation  were  taken  to  be  the  steady-state  values  obtained 
via  the  isothermal  model.  These  variables  will  not  satisfy  the  equations  of  the  more  general 
model  initially,  and  will  evolve  in  time  until  a  new  steady  solution  is  reached.  A  comparison 
of  the  two  steady  solutions  should  then  provide  an  indication  of  how  accurate  the  isothermal 
model  is. 

For  the  isothermal  case,  we  showed  in  the  previous  section  that  no  steady-state  solution 
existed  for  w{z  =  ()  greater  than  C,o  where  C,o  rcpiescnted  the  isothermal  sound  speed.  For 
u’(r  =  f)  —  C,o,  the  flow  is  choked  and  it  is  easy  to  show  that  the  breech  and  exit  pressures 
must  satisfy 


To  =  P(2=0)=2P(2  =  f).  (19) 

Under  conditions  which  are  not  isothermal,  the  sound  speed  is  more  difficult  to  compute,  and 
we  have  for  simplicity  used  Eq.  (49)  to  supply  the  additional  boundary  condition  required. 
The  sound  speed  at  the  point  z  =  C  has  been  estimated  in  our  calculations  from  the  very 
approximate  expression 


Cs  =  (50) 

where  ')efj  is  the  effective  adiabatic  exponent  or  specific-heat  ratio  for  the  plcisma  given  by 
(Zeldovich  and  Raizer  19CG) 


7c//  =  1  -f  PlipE). 


(51) 


Equation  (50)  is  known  (Tidman  and  Goldstein  1985)  to  produce  an  estimate  of  Cs  which 
is  slightly  higher  than  that  produced  by  the  more  rigorous  expression,  namely, 


-'S), 


1/2 


(52) 


with  the  subscript  S  denoting  partial  differentiation  at  constant  entropy. 
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Silowii  in  F'igs.  11  and  12  arc?  pressure  and  temperature  profiles  obtained  from  bolli  the 
isotlicrnial  model  and  the  more  general  model  discussed  above.  Hermits  are  plotted  as  a 
function  of  the  dimensionless  distance  ^  =  zjC.  As  can  be  seen,  the  pressure  in  the  arc 
agrees  to  within,  at  worst,  about  1 5 everywhere;  the  temperature  is  in  somewhat  better 
agreement.  Presumably,  the  slight  drojj  in  lernperature  with  increasing  if  results  from  a 
greater  |)crccntagc  of  the  total  energy  in  the  arc  being  as.sociated  with  kinetic  energy  of  the 
gas  products  as  they  are  accelerated  near  exit.  Similar  observations  were  noted  in  PZl.  The 
load  voltages  predicted  by  the  general  and  isothermal  (see  Table  6)  models  were  20‘17  and 
21SS  respectively.  Similarly  reasonable  agreement  has  been  found  for  other  variables  and 
other  shots  and  we  conclude  that  the  isothermal  model  can  be  used  reliably  for  predictions 
in  cases  for  which  the  flow  is  steady  and  for  which  e.xtremc  accuracy  is  not  required. 


0  0.2  C.4  0.6  C.8  1 

Figure  11.  Comparison  of  Steady-State  Pressure  Profiles  in  the  Isothermal  and  General 
Models  for  Shot  8. 

6.2  .  Time-Dependent  Calculations.  In  addition  to  the  steady-state  calculations  de¬ 
scribed  above,  we  have  also  undertaken  a  number  of  time-dependent  calculations  with  the 
more  general  model.  The  purpose  of  these  calculations  was  to  compare  with  time-dependent 
experimental  data  for  the  shots  discussed  in  the  preceding  subsection.  It  is  necessary  for 
these  calculations  to  input  the  current  j(f)  as  a  function  of  time.  This  current  was  obtained 
from  experiment  (see,  e.g.,  Fig.  5)  and  then  fitted  to  a  sixth-order  polynomial  in  tirrie.  Such 
a  fit  was  found  to  give  an  accurate  representation  of  the  actual  measured  current  except  at 
the  very  end  and  the  very  beginning  of  the  current  pulse.  These  extremes  were  omitted  from 


36 


0  e  r  c " I  3 1  j  c  .  S  ?  z  •  g 


^  r 


C.2  0.4  0.6  C.5  1 

Figure  12.  Comparison  of  Steady-State  Temperature  Profiles  in  the  Isothermal  and  General 
Models  for  Shot  8. 

the  calculation  anyway,  since  we  selected  for  analysis  a  region  about  peak  current  for  which 
the  pressure  measurements  were  felt  to  be  particularly  accurate. 

For  all  these  calculations,  the  boundary  condition  at  exit  was  taken  to  be  that  given  in 
Eq.  (49).  When  this  condition  was  applied,  we  generally  found  that  the  exit  speed  of  the 
gas  products  was  slightly  less  than  our  estimate  of  the  sound  speed  at  exit  obtained  from 
Eq.  (50).  The  initial  conditions  employed  in  each  case  were  the  steady-state  values  of  the 
various  parameters  obtained  at  a  value  of  the  current  equal  to  that  which  we  chose  to  be  the 
initial  value  in  the  pulse  described  above.  These  values  were  obtained  by  running  the  lime- 
dependent  code  at  the  initial  current  until  a  steady  stale  was  achieved,  and  then  allowing 
the  current  to  vary. 

We  now  discuss  some  of  the  results  of  the  calculations.  We  will  confine  our  attention 
to  Shots  4  and  8;  other  cases  gave  very  similar  results.  We  will  first  discuss  some  purely 
theoretical  results  and  then  compare  with  the  available  experimental  data. 

Shown  in  Figs.  13-15  are  pressure,  velocity,  and  temperature  for  Shot  8  at  three  different 
times.  The  only  noteworthy  eispect  of  these  graphs  is  the  observation  that  the  profiles  look 
reasonably  steady;  the  shapes  of  the  curves  hardly  change  as  the  current  rises  from  about 
3G  kA  at  300  fis  to  a  peak  of  about  58  kA  at  around  460  /rs,  and  then  drops  to  about  26  kA 
at  700  fis.  Further  evidence  of  the  steady  behavior  can  be  seen  from  data  in  Table  7  in  which 
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are  shown  in  Columns  3  and  4  the  power  dissipated  in  th.c  arc  at  tlu’  time  in  cjuestion.  Vo- 
and  the  power  passing  tlie  capillar}'  exit  plane.  Vy.  'I'lie  latter  (luantit}'  wa.-'  calculated  from 
the  expression 


Vy  =  A[pw(E  -r  W'/2)  +  p3) 

As  suggested  previously,  one  would  expect  results  in  Columns  3  and  4  to  be  equal  in  steady 
state.  Actually  the  results  agree  to  within,  at  worst,  about  5%  .*  .Also  shown  in  the  table 
in  Column  6  is  the  estimated  sound  speed  at  the  exit  plane  calculated  from  Eq.  (50)  and, 
in  Column  6,  the  exit  velocity  of  the  gaseous  products.  .As  pointed  out,  the  exit  speeds  are 
somewhat  smaller  than  the  estimated  sound  speed. 
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Figure  13.  Calculated  Pressure  Profiles  at  Several  Times  for  Shot  8. 

One  of  the  principal  reasons  for  investigating  a  time-dependent  model  was  the  obser¬ 
vation,  made  by  us  as  well  as  others,  that  many  of  the  plasmadynamic  variables  in  these 
arc-discharge  experiments  exhibit  significant  time-dependent  behavior.  The  very  early  ex¬ 
periments  of  Ibrahim  (1980),  for  example,  indicated  that  the  capillary  voltage  was  not  a 
single-valued  function  of  the  current  as  one  would  expect  for  completely  steady  behavior. 
Furthermore,  it  Wcis  noted  in  that  work  that  even  more  significant  time-dependent  behavior 
was  observed  in  other  arc  variables  such  &s  the  pressure.  The  later  work  by  Ruchti  and 

*We  could  also  comp£U'e  the  mass  ablated  per  unit  time  with  the  mass  which  flows  past  the  exit,  as  well  the  momentum 
fljx  to  the  total  pressure  drop  across  the  capilltuy.  These  calculations  have  in  fact  been  carried  out,  and  the  results  indicate 
nearly  steady  behavior. 
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Table  7.  Theoretical  Results  from  Time-Dependent  Calculations  for  Shot  S. 


Cs{(=l) 

(^•77l/.<.) 

r(,^=  1) 

(  A'77|/s) 

300 

36.1 

60.7 

9.8 

8.S 

500 

58.3 

121.2 

121.7 

11.1 

10.0 

700 

26.2 

38.5 

41.2 

8.9 

8.2 

Niemeyer  (1978),  if  examined  carefully,  also  demonstrates  substantial  time  dependence  in 
the  measured  pressure  versus  current,  with  the  current  reaching  peak  value  considerably 
prior  to  the  time  at  which  the  pressure  peaks.  Finally,  with  regard  to  electrothermal  guns, 
"hysteresis”  effects  have  been  observed  in  voltage-versus-currenl  profiles  not  only  by  us  but 
by  other  investigators  at  BRL  as  well  (Katulka  ei  al.  1990). 

In  order  to  compare  calculated  time-dependent  results  with  experimental  data,  we  have 
plotted  in  Figs.  16-19  for  both  Shots  4  and  8  the  experimentally  measured  conductance  and 
pressure  as  a  function  of  current.  The  graphs  were  plotted  only  for  the  time  during  which  the 
calculations  were  carried  out,  i.  e.,  for  200  fis  <  t  <  700 /is  for  Shot  8,  and  300  fis  <  t  <  700/is 
for  Shot  4.  Shown  on  the  same  graphs  are  the  calculated  results  for  the  appropriate  time 
interval. 

In  Fig.  16  are  the  experimental  and  theoretical  values  of  the  conductance  for  Shot  8  as  a 
function  of  current.  It  is  evident  from  the  curves  that  both  the  experimental  and  theoretical 
value  of  the  conductance  is,  at  any  fixed  value  of  i,  slightly  lower  when  the  current  is 
rising  than  when  it  is  falling.  The  difference  is  much  more  pronounced,  however,  in  the 
experimental  than  in  the  theoretical  case.  Similarly,  shown  in  Fig.  IT,  are  the  calculated  and 
experimental  values  of  the  pressure  Pp  at  the  point  ^  =  0.28  (or  z  =  17  mm)  as  a  function 
of  the  input  current.  Here  the  time-dependent  effect  is  considerably  more  evident,  especially 
in  the  experimental  data,  where  the  pressure  is  a  factor  of  two  or  so  higher  when  i  is  falling 
compared  to  when  it  is  rising.  Similar  results  can  be  seen  for  Shot  4  in  Figs.  18  and  19. 

The  slight  time  dependence  in  the  theoretical  results,  which  arises  from  the  time- varying 
current,  can  be  described  approximately  as  follows.  A  small  increase  in  current  leads  to  an 
immediate  increase  in  Joule  heating  and  a  corresponding  rise  in  temperature.  Thus,  the 
rate  of  ablation,  which  depends  primarily  on  T,  also  rises  with  the  temperature.  The  rate 
of  mass  loss  in  the  capillary,  however,  depends  not  only  on  T  but  also  varies  strongly  with 
the  density  p.  Therefore,  the  rate  at  which  mass  flows  out  the  capillary  will  be  lower  than 
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Figure  16.  Experimental  and  Calculated  Values  of  the  Conduct<ince  as  a  Function  of 
Current  for  Shot  8. 
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Figure  17.  ExperinientaJ  and  Calculated  Values  of  the  Pressure  as  a  Function  of  Current 
for  Shot  8. 
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Figure  18,  Experimental  and  Calculated  Values  of  the  Conductance  a«  a  Function  of 
Current  for  Shot  4. 


•  Shol  4  E*p*rimenl  --  Cufftnl  riling 
®  Shot  4  E«p#rim*nf  --  Currir^l  lolling 
Shot  4  Thiofy  —  Currint  filing 
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Figure  19.  Experimental  and  Calculated  Values  of  the  Pressure  as  a  Function  of  Current 
for  Shot  4. 
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tlif  al  wliltli  ll  llicffn^cs  owiii^  lo  nliinliiHi,  aiul  Up! Hi  will  lie  ivc,  Ah  A  r»'hiill, 
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which  liil|hl  llliimlliAle  0111  niidei alalnlilig. 

O.fl  i 'ftielllfliloUl  Ld.liploylim  Nlyrtainytl  Preaaiiie  na  liijnit.  One  poaaihllllv  for 
ill)  lack  of  Agieciiienl  helMccli  (In  eapetliiienlal  ainl  iheorellial  salne*  of  llie  pieaalire  c'uiihl 
hi  lliai  (In  tall  ipud  of  (In  gniet  la  ool  Well  np|noalmAli')|  hy  (In'  aoiih  condition.  In 
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llcsults  of  the  calculation  for  Sliot  8  arc  shown  in  Fig.  '20  in  which  arc  i)lotlcd  the 
plasma  conductance  oblaincil  from  the  experimentally  incnsurccl  voltage  and  current,  the 
value  obtained  via  the  steady-state  model  under  the  as.sumption  that  the  flow  was  sonic, 
and  finally  the  value  obtained  when  the  calculation  w;is  carried  out  using  the  experimental 
pressure  as  descrihed  above.  'I  he  results  tend  to  indicate  slightly  heller  agreement  between 
theory  and  experiment  for  the  case  in  which  the  pre.sHure  was  matched  than  in  the  case 
in  which  the  exit  flow  was  assumed  to  be  sonic.  As  is  evident  in  the  figure,  and  as  has 
hern  pointed  out  previously,  the  voltage  or  conductance  is  not  very  sensitive  to  the  capillary 
pressure  since  the  jjlasrna  conductivity  depends  primarily  on  the  tcmi)erature, 


Figure  20,  Measured  and  Tlieorelical  Values  of  the  Oondurlance  for  Sonic  and  Nonsonic 
I'lxit-l'luw  Conditions,  'I'he  results  are  for  Shot  ft, 

It  Is  also  of  Interest  to  determine  the  extent  to  which  the  sonic  coiulition  is  satisfied 
when  the  plesBUie  is  mat  died  t  o  experiment  in  the  manner  d(  serihed  above,  'I'o  answer  this 
qiii'slion  we  have  plotted  the  .Mach  miinher  of  the  gas  at  exit  for  the  lime  during  which 
the  falciilation  was  carried  out.  Hesulls  are  shown  in  Fig,  21,  As  is  evident  the  flow  is 
always  luhsonic  since  the  Mach  rniiiilM'r  is  less  than  unity.  Ibr  the  first  d,*)!)  /i.s  or  so  of 
the  <  alciilatioii,  the  gas  exits  the  tube  at  a  roughly  eoiistanl  fraction  of  the  sound  speed  . 
AroiJiiil  the  lime  of  current  peah,  however,  the  Math  iiuinher  begins  to  decrease  rapidly  and 
dioj»s  by  a  factoi  «>f  about  two  during  the  remainder  of  the  t  ule'jtalioii.  These  results  will 
he  tlisniswed  ill  some  giealer  detail  in  the  following  set  lion. 
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Figure  21.  Plasma  Mach  Number  at  Exit  Plane  for  the  Pressure-Match  Calculation. 

7.  CONCLUSIONS  AND  DISCUSSION 

We  have  carried  out  both  steady-state  and  some  time-dependent  calculations  of  plasma 
properties  for  capillary  discharges  similar  to  those  which  occur  in  electrothermal  guns.  Re¬ 
sults  have  been  compared  with  experimental  values  of  the  capillary  voltage  and  of  the  pres¬ 
sure  measured  at  a  specific  point  in  the  capillary  tube.  General  conclusions  obtained  from 
this  work  can  be  described  as  follows: 

1.  The  steady-state,  isothermal  model,  for  the  assumption  that  the  flow  is  sonic  at  exit,  gives 
estimates  of  the  plasma  voltage  which  arc  in  reasonable  agreement  with  the  experimental 
values,  The  voltage  appears  to  be  fairly  steady  although  some  hysteresis  is  evident  in  the 
experimental  voltagc-versus-current  curves, 

2.  The  measured  and  theoretical  values  of  the  pressure  do  not  agree,  differing  typically  by 
as  much  as  a  factor  of  two.  Substantial  time-dependent  behavior  of  the  pressure  is  indi¬ 
cated  from  experimental  prcssure-versus-currcnt  traces.  The  time-dependent  model,  which 
accounts  for  the  time  variation  of  the  current  profile,  predicts  the  correct  trends  in  the  non¬ 
stationary  behavior  of  both  the  voltage  and  pressure,  but  the  magnitude  of  the  effects  is 
significantly  underestimated. 

3.  Calculations  in  which  the  flow  is  not  assumed  to  be  sonic  at  exit  but,  rather,  the  exit 
speed  is  determined  by  using  the  exi)erimental  value  of  the  pressure,  were  also  carried  out. 
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Results  of  these  calculations  indicate  slightly,  but  unconvincingly,  better  agreement  with 
experimental  voltage  meeisurements  than  do  calculations  performed  for  the  sonic  condition. 

Although  we  are  reduced  largely  to  speculation,  we  will  now  discuss  a  number  of  reasons 
which  might  contribute  to  the  discrepancy  between  theory  and  experiment  <md  suggest  future 
improvements  in  both  the  models  and  experiments  that  could  illuminate  the  subject  further. 

The  first  and  most  obvious  limitation  of  using  the  model  to  compare  with  experimental 
results  is  the  necessity  of  assuming  some  boundary  condition  concerning  the  exit  velocity  of 
the  gaseous  products.  Fortunately,  many  of  the  plasma  properties,  such  as  its  conductivity, 
do  not  depend  strongly  on  this  condition.  One  possibility  for  the  nonsteady  behavior  of  the 
capillary  properties  could  be  a  rise  in  pressure  in  the  brass  extender  tube,  resulting  in  an 
unchoking  of  the  flow  as  time  progresses.  Evidently,  however,  no  such  effect  could  account 
for  the  time-dependent  voltage  observed  when  capillaries  are  fired  into  open  air,  such  as 
in  the  experiments  done  by  Katulka,  Burden,  and  White  (1990),  and  other  explanations 
must  be  investigated.  One  conceivable  explanation  could  be  the  substantial  amount  of 
electrode  ablation  observed  near  the  exit  of  capillary.  This  effect  might  be  expected  to 
lead  to  a  situation  in  which  the  heavy,  less  mobile,  electrode  material  prevents  the  escape 
of  the  capillary  plasma.  If  such  an  effect  does  occur,  it  would  become  more  dominant  with 
increasing  time  since  the  ablation  of  metallic  material  will  lag  the  ablation  of  polyethylene  by 
a  time  period  which  depends  very  strongly  on  the  temperature.  This  observation  is  consistent 
with  the  results.  Electrode  ablation  could  be  included  in  the  models  discussed  above,  but  the 
analysis  would  require  the  treatment  of  a  non  homogeneous  mixture  of  polyethylene  and  the 
electrode  material.  We  are  currently  attempting  to  develop  such  a  model.  Experimentally, 
the  escape  velocity  of  the  gas  at  capillary  exit  could  be  measured.  We  have  in  fact  attempted 
these  measurements  via  use  of  fiber-optic  cables  placed  close  together  near  the  tube  exit. 
Metisurements  similar  to  these  have  been  performed  successful!}-  in  the  past  in  the  analysis 
of  railgun  plasma  armatures  (Zielinski  and  Powell  1990).  For  the  cases  investigated  here, 
however,  and  particularly  for  the  higher-current  shots,  the  large  magnitudes  of  the  velocities 
involved  have  made  the  data  difficult  to  evaluate.  Much  of  the  difficulty  arises  because  the 
technique  provides  only  distance-versus-time  data  and  numerical  differentiation,  subject  to 
significant  error,  is  necessary  to  obtain  the  velocity.  Preliminary  results  have  suggested, 
however,  that  the  exit  velocities  are  considerably  less  than  the  local  sound  speed  in  the 
gas  at  the  point  in  question.  Further  experimental  work  is  required,  however,  before  any 
quantitative  statements  can  be  made.  It  would  also  be  worthwhile  to  undertake  experiments 
with  different  types  of  electrode  materials  in  an  effort  to  determine  the  extent  to  which  the 
electrodes  affect  the  results. 

A  second  reason  for  the  lack  of  agreement  between  theory  and  experiment  is  the  possi- 
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bility  that  the  capillary  plasma  may  be  ’’contaminated.”  As  indicated  previously,  this  and 
other  capillary-discharge  models  are  based  on  the  assumption  that  the  capillary  consists  of  a 
mixture  of  the  various  species  which  comprise  the  capillary  wall.  Clearly,  however,  if  there  is 
significant  electrode  erosion,  whether  at  the  breech  or  muzzle  end  of  the  discharge  tube,  the 
electrode  material  will  become  a  part  of  the  entrained  plasma.  The  resulting  mixture  may 
then  have  properties  different  from  those  of  the  assumed  mixture.  Furthermore,  particularly 
if  the  flow  is  not  choked,  there  may  be  species  diffusion  from  the  region  beyond  the  breech 
back  into  the  capillary.  Again,  under  this  condition,  the  constituency  assumed  for  the  plasma 
will  not  be  correct  and  its  properties  will  be  modified  in  some  unknown  way.  We  expect  that 
this  contamination  may  affect  significantly  sensitive  plasma  variables,  such  as  the  pressure, 
density,  and  velocity,  but  not  insensitive  ones  such  as  the  voltage  or  temperature. 

A  third  possibility  for  the  lack  of  agreement  may  be  that  the  current  is  not  being  con¬ 
ducted  throughout  the  entire  radius  of  the  capillary  or,  in  other  words,  the  thickness  of  the 
vapor  zone  may  not  be  negligible  as  assumed.  In  his  early  experiments,  Ibrahim  (1980)  found 
it  necessary  to  postulate  such  an  effect  in  order  to  achieve  an  energy  balance  between  the 
radiated  and  dissipated  power.  Perhaps  these  or  similar  conclusions  have  been  responsible 
for  the  development  of  a  number  of  "two-zone”  models  such  as  proposed  by  Kovitya  and 
Lowke  (1984).  In  tho.se  models  the  vapor  layer  is  treated  as  a  region  of  finite  thickness 
whose  dimensions  vary  with  the  properties  of  the  plasma  in  the  central  core.  Necessarily,  the 
treatment  of  the  vapor  layer  is  highly  approximate  and  a  number  of  assumptions,  difficult  to 
justify,  are  made  concerning  its  properties.  Nevertheless,  the  models  discussed  in  this  report 
could  presumably  be  extended  to  two  zones.  It  is  evident  that  conduction  of  the  current  in  a 
region  having  cross  sectional  area  smaller  than  the  tube  itself,  will  produce  higher  pressures 
than  for  the  case  when  the  conduction  occurs  across  the  entire  cross  section.  This  aspect  of 
the  results  would  clearly  be  in  better  agreement  with  the  experimental  data.  However,  the 
two- zone  models  referred  to  above  have  all  been  steady  state  and  are  clearly  not  capable  of 
explaining  the  time-dependent  behavior  of  the  pressure-  or  voltage-versus-time  curves.  An 
extension  of  those  models  to  the  time-dependent  case  would  be  necessary  before  one  could 
hope  to  obtain  any  explanation  of  this  aspect  of  the  problem.  Of  course,  the  two-zone  models 
are  just  a  rough  approximation  to  a  real  two-dimensional  model.  As  pointed  out  in  PZl, 
a  time-dependent,  two-dimensional  investigation  of  this  problem  is  appropriate  and  timely. 
Experimentally,  the  radius  of  the  current-carrying  part  of  the  capillary  could  probably  be 
inferred  from  measurements  of  the  magnetic-induction  field  as  a  function  of  radial  distance 
in  the  tube.  Somewhat  similar  experiments  (Zielinski  1989b)  have  been  performed  in  the 
past  in  railguns. 

A  fourth  reason  could  be  the  experiment  itself.  Although  the  pressure  transducer  was 
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calibrated  in  situ,  the  time  scale  of  the  calibration  was  an  order  of  magnitude  greater  than 
that  of  the  experiment.  Ordinarily,  this  difference  would  not  be  significant,  but  since  a 
polyethylene  wall  shields  the  gage,  the  dynamic  properties  of  the  polyethylene  may  affect 
the  results.  In  addition,  it  would  be  helpful  to  measure  the  pressure  at  a  number  of  points 
along  the  length  of  the  capillary,  and  for  capillaries  of  fixed  diameter. 

A  final  reeison  that  theory  and  experiment  may  not  agree  completely  is  that  there  may 
exist  mechanisms,  other  than  the  thermal  one  assumed,  which  contribute  to  removal  of  wall 
material.  Such  mechanisms  include,  for  example,  pyrolysis  and  stress-wave  effects  which  can 
cause  spallation.  These  mechanisms,  and  the  range  of  incident  power  densities  for  which 
they  are  important,  have  been  discussed  by  Nemes  and  Randles  (1989)  and  are  important 
for  these  experiments. 
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